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STOCHASTIC PARTICLE APPROXIMATION OF THE KELLER-SEGEL 
EQUATION AND TWO-DIMENSIONAL GENERALIZATION OF BESSEL 

PROGESSES 

NICOLAS FOURNIER AND BENJAMIN JOURDAIN 


Abstract. The Keller-Segel partial differential equation is a two-dimensional model for chemo- 
taxis. When the total mass of the initial density is one, it is known to exhibit blow-up in finite 
time as soon as the sensitivity x of bacteria to the chemo-attractant is larger than Stt. We 
investigate its approximation by a system of N two-dimensional Brownian particles interacting 
through a singular attractive kernel in the drift term. 

In the very subcritical case x < ^tt, the diffusion strongly dominates this singular drift: we 
obtain existence for the particle system and prove that its flow of empirical measures converges, 
a,s N oo and up to extraction of a subsequence, to a weak solution of the Keller-Segel equation. 

We also show that for any N > 2 and any value of x > 0? pairs of particles do collide with 
positive probability: the singularity of the drift is indeed visited. Nevertheless, when x < 27rN, 
it is possible to control the drift and obtain existence of the particle system until the first time 
when at least three particles collide. We check that this time is a.s. infinite, so that global 
existence holds for the particle system, if and only if x ^ 87r(A’ — 2)/(N — 1). 

Finally, we remark that in the system with N = 2 particles, the difference between the two 
positions provides a natural two-dimensional generalization of Bessel processes, which we study 
in details. 


1. Introduction and results 


1.1. The model. The Keller-Segel equation, introduced by Patlak [30] and Keller and Segel [2T| . 
is a model for chemotaxis. It describes the collective motion of cells which are attracted by a 
chemical substance and are able to emit it. In its simplest form it is a conservative drift/diffusion 
equation for the density ft{x) > 0 of cells (particles) with position x G at time t > 0 coupled 
with an elliptic equation for the chemo-attractant concentration. By making the chemo-attractant 
concentration explicit in terms of the cell density, one obtains the following closed equation; 

(1) dtft{x) + xdivxiiK * ft){x)ft{x)) = A^ft{x), 

where x > 0 is the sensitivity of cells to the chemo-attractant and where 


( 2 ) 


K{x) = 


27r|xp 

In the whole paper, we adopt the convention that K{0) = 0. 

This equation preserves mass and ft{x)/J^ 2 foiy)dy solves the same equation with x replaced 
by xJg 2 fo{y)dy. We thus may assume without loss of generality that J^ 2 foix)dx = 1. 
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As is well-known, we have formally ^ xft{x)dx = 0 and ^ \x\^ ft{x)dx = 4 — x/(27r). 

Consequently, introducing Vt := /r 2 |a; - /^a yft[y)dy\^ft{x)dx, it holds that -^Vt = 4 - x/(27r). 
Since Vj is nonnegative, some kind of blow-up necessarily occurs before time 27rVo /(x — Stt) when 
X is larger than the critical value Stt. 

Concerning the well-posedness theory, let us mention Jager and Luckhaus m, Blanchet, Dol- 
beault and Perthame [1], Dolbeault and Schmeiser [4] and Egaha and Mischler [5]. In particular, 
the existence of solutions is verified in [16] (for sufficiently smooth initial conditions), these solu¬ 
tions being local (in time) if x > 0 is large and global if x > 0 is small. The existence of a unique 
strong (in some precise sense) solution when x < Stt is shown in [I] (existence) and [^ (uniqueness), 
still for reasonable initial conditions. The main tool is the free energy and its relation with its time 
derivative. By passing to the limit in a sequence of regularized Keller-Segel equations where the 
kernel K is replaced by a bounded kernel and by introducing defect measures to take into account 
blow-up, the existence of generalized weak solutions to 0 is checked in |4], even when x > Stt. 
The blow-up phenomenon has been investigated by Herrero and Velazquez [niMiinz]- We refer 
to Horstmann mils] and Perthame m for review papers on this model. 

1.2. Weak solutions. We denote by 7^(IR^) the set of probability measures on and we set 
7^1 (K^) = {/ £ P(R^) : TOi(/) < oo}, where mi(/) = J^ 2 \x\f{dx). We will use the following 
notion of weak solutions. 


Definition 1. Let x > 0 T £ (0,oo] be fixed. We say that a measurable family {ft)t^io,T) of 
probability measures on is a weak solution to o on [0,T) if the following conditions hold true: 

(a) for all t £ [0,T), J* - y\-^fsidy)fsidx)ds < oo; 

(b) for all (j) £ C2(R2)^ all t £ [0,r), 


/ 

Jr: 


(l){x)ftidx) = [ (j){x)fo(dx) + [ [ A(j){x)fs{dx)ds 

JR2 Jo JR2 

+ X [ f [ K{x - y) ■V4){x)fs{dy)fs{dx)ds. 
Jo JR2JR2 


Of course, (a) implies that everything makes sense in (b). Performing a symmetrization in 
the last term leads to another weak formulation of 0 which requires less stringent integrability 
conditions, but which is not suitable in view of the following probabilistic interpretation. 

1.3. The associated trajectories. We now introduce a natural probabilistic interpretation of 
the Keller-Segel equation. 

Definition 2. Let x > 0 and T £ (0,oo] be fixed. We say that a M.'^-valued continuous process 
iXt) te[o.r) adapted to some filtration {J^t)t^[o,T) solves the nonlinear SDE ([Sj) on [0,r) if, for 
ft ■= c(Xt), it holds that 

(°-) /o/R 2 /R 2 |a: - y\~'^f,idy)fsidx)ds < oo for all t £ [0,T); 

(b) there is a 2-dimensional {!Ft)t^[o,T)~Bxownian motion (St)tg[o,T) such that for all t £ [0,T) 


(3) 


Xt = Xo + V2Bt 


X [ {K * fs){X,)ds. 
Jo 


The main idea is that (Vt)tg[o^T) represents the time-evolution of the position of a typical cell, 
in an infinite system of cells undergoing the dynamics prescribed by the Keller-Segel equation. The 
following remark immediately follows from the ltd formula. 
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Remark 3. Let x > 0 6e fixed. For (-^t)tG[o,T) solving the nonlinear SDE Q, the family (/* = 
F{Xt))t^\Q^T) 'Is 0 . weak solution to the Keller-Segel equation (fT|). 

1.4. The particle system. We next consider a natural discretization of the nonlinear SDE: we 
consider N > 2 particles (cells) with positions ,... ,X^'^ solving (recall that A"(0) = 0) 

V ^ 

(4) Xf^ = + + K{Xl’^ - Xfi^)ds. 

i=i 

More precisely, a solution on [0,r) is a continuous (M^)^-valued process 7 v,tG[o,T) 

adapted to some filtration (J^t)tg[o^T) if the initial conditions Xq, * = 1,... ,iV are i.i.d. with com¬ 
mon law /o € P(M^) and if there is a 2Wdimensional (J^i)tg[o_T)-Brownian motion ..., B^)t>o 

such that O holds true for all t G [0, T) and alH = 1,..., A^. 

Of course, such a particle system is not clearly well-defined, due to the singularity of K. More¬ 
over, the singularity is visited, as shown by the following statement. 

Proposition 4. For any N > 2, any x > 0, any fo G P(R^), any to > 0 and any solution (if it 
exists) to dll), 

P(3 s G [0,to], ^l<i<j<N : X)’^ = > 0. 

However, we expect that particles are almost independent (for N large) and look like N copies 
of the solution to the nonlinear SDE, at least in the subcritical case x S (0, Stt) or locally in time 
in the supercritical case x ^ Stt. This problem seems important, both from a physical point of 
view, as a step to the rigorous derivation of the Keller-Segel equation, and from a numerical point 
of view. 


1.5. Main resnlts. We first check that the particle system ((ll) exists when x is (very) subcritical. 

Theorem 5. Let N >2 and x G (0, 2TrN/{N — 1)) be fixed, as well as fo G Pi(R^). There exists 
a solution {Xl’^)t^[o^ao),i=i,...,N to (jll). Furthermore, the family {(X(’^)tg[o_oo)! * = is 

exchangeable and for any a G {{N — 1)x/{‘2ttN), 1), any T > 0, 

i-T 


(5) 


E 




2,N\a-2 


ds 


< 


{2V2{fo,^l + \x\-^)+4V2Ty 

a(2a -{N - l)x/(7riV)) 


As already mentioned, such a result is not obvious, since K is singular and since its singularity 
is visited. The main point is to observe that dS]) a priori holds true for some a < 1. This will 
imply that that E[|Ar(X]^’^ — Xg’^)|] should be controlled (with some margin since a — 2 < —1). 
This will be sufficient to prove existence by compactness. The formal computation is as follows: 
by the Ito formula, for a G (0,1), 

= V2a\Xl'^-X^ 


( 6 ) 


d\Xl'^-Xf'^\ 


-Xy^y-{dBl-dBt) 


+ 2a'^\xl’^ - xy^y-^dt -- xy^y-^dt 

ttN 


ax 

2ttN 


\xl’^ -xy^y-\xl'^ -X, 


2,N- 


N 

-n 

=3 


xr-x, 


1,N 


X, 


2,N 


-XI 


,Af 




.1,A| 


\XfX _ X)’ 


N,' 


]dt. 


The second term in the right-hand side is the Ito correction due to diffusion, the third term is the 
contribution of the interaction between the particles 1 and 2 and the last term is the contribution 
of the interactions with between particles 1,2 and the rest of the system. By exchangeability 
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and Holder’s inequality, the expectation of the last term in the right-hand side is greater than 
— [a(-/V — 2)x/(7r7V)]E[|X(^’'^ — \ The assumption y < 2ttN/{N — 1) ensures us that the 

Ito correction dominates the drift contribution. More precisely choosing a G (x(iV— l)/(27rA^), 1), 
integrating in time and taking expectations, one obtains 

a(2a - - Xyr-^]ds < E[\Xl’^ - Xf’^n 

The right-hand side is easily bounded, uniformly in N, using the oddness of K, whence (l5|). A 
similar computation was performed by Osada in |291 Lemma 3.2] for systems of stochastic vortices. 

Next, and this is the main result of the paper, we show some tightness/consistency as —>■ oo in 
the (very) subcritical case x < 27r. Such a result follows quite easily from the the bound ([S]), which 
is uniform in N (when x < 27r). We endow C([0, oo), R^) with the topology of uniform convergence 
on compact time intervals, and ^((^([O, oo), R^)) with the associated weak convergence topology. 
Finally, we endow (^([O, oo), P(R^)) with the topology of uniform convergence on compact time 
intervals associated with the weak convergence topology in P(R^). 

Theorem 6. Let x G (0, 27r) be fixed, as well as fo G For each N > 2, consider the 

particle system (Aj’^)tg[Q oo),i=i,...,Ar built in Theorem\^ as well as the empirical measure = 

which a.s. belongs to P(C'([0,oo),R^)). For each t>0, we also set = 

which a.s. belongs toV(JSfi). 

(i) The sequence {/i^, N > 2} is tight m ^((^([O, oo), R^)). 

(ii) Any (possibly random) weak limit point p, of {h^)n >2 is a.s. the law of a solution to the 
nonlinear SDE 0 with initial law /q. 

(Hi) In particular, we can find a subsequence such that (pL^'‘)t>o goes in law, as k ^ oo, in 
(^([O, oo),P(R^)), to some {nt)t>o, which is a.s. a weak solution to o starting from p-o = fo. 

We are quite satisfied, since this result seems to be the first result concerning the convergence of 
the true particle system (without cutoff) to the Keller-Segel equation. However, there are two main 
limitations. First, this result should more or less always hold true in the subcritical case x G (0, Stt). 
Second, we are not able to prove the convergence, we have only compactness/consistency. This is 
due to the fact that we are not able to prove that our limit point (/it)t>o a.s. belongs to the class 
of weak solutions in which uniqueness is known to hold true. Thanks to Egaha and Mischler 0 , it 
would suffice to show that (/it)t>o satisfies the free energy dissipation inequality, which is slightly 
stronger than the requirement (/rt)t>o S np>iL;^gj,([0,oo),LP(R^)) a.s. We believe this is a very 
difficult problem. 

We next prove that, when x < ^ttN, the particle system always exists until 3 particles encounter. 
In view of 0, this is not surprising. Indeed, the assumption x < SttA ensures us that the 
Ito correction still dominates the contribution of the interaction between the particles 1 and 2. 
Moreover, it is not very hard to control the last term of 0 until a 3-particle collision occurs. 

Theorem 7. Let x > 0, A > max{2, x/(27r)} be fixed, as well as fo G Pi(R^) such that fo{{x}) = 
0 for all X G R^. There exists a solution (Aj’^)tg[o_T;v).i=i,....Ar to (0), with 

ttv = sup inf 41 > 0 : 3 i,j,k pairwise different such that 

i>i '■ 

\X^,N _ ^ ^ 
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The family * = 1, • ■ •, N} is exchangeable and for any a G (x/(27rA^), 1), 

(7) a.s., for all t G [0,tn), f \Xl’^ — Xg’^\°‘~'^ds < oo. 

Jo 

Finally, it holds that (i) tn = oo a.s. if x ^ 8'!r{N — 2)/(7V — 1) and (ii) Tpf < oo a.s. if 
X>87r(iV-2)/(7V-l). 

This result thus in particular shows the global existence for the particle system in the subcritical 
case X < Stt for all N large enough. This result seems to be new, as well as our method to check 
it, which is quite specific to the model. 

As we will see in the proof of Lemma ITSl-SteD 2, for any subsystem I C {l,...,Af}, the process 
Rj. = ~ where Xf = \I\~^ ^ behaves like the square of a Bessel 

process of dimension (|/| —1)(2 — (x|7|)/(47rAf)), when neglecting the contribution of the interaction 
with the other particles. Similar computations for I = {!,...,A^} were performed by Haskovec 
and Schmeiser [TUI Page 139] and Fatkullin [UJ Page 89]. The condition x < — 2)/(N — 1) 

implies that for all |/| = 3,..., N, the dimension (|/| — 1)(2 — (x|7|)/(47rA^)) is greater than 2, so 
that does never reach 0: there are no collisions involving more than two particles. Of course, 
the situation is actually much more complicated, since we have to justify that of all I, we can 
indeed neglect the contribution of the interaction with the other particles. 

Remark 8. When x G (0,87r), we thus show that, for N large enough, the particles labelled 
1,2,3 do a.s. never encounter. To extend the tightness/consistency result of Theorem\^to some 
X G [27r,87r), we believe that a quantitative and uniform (in N) version of this fact might be 
sufficient. 

Finally, we study the case of two particles N = 2. The average of the two positions is a two- 
dimensional Brownian motion and their difference Dt follows an autonomous SDE with singular 
drift driven by a Brownian motion, which can be seen as a natural two-dimensional generalization 
of a Bessel process of dimension (2 — x/(47r)). We show that the equation for Dt (as well as (|T|)) is 
nonsense when x > 4?!, in that there cannot exist global solutions. But this is only a small problem 
related to the fact that a Bessel process with dimension 5 G (0,1] does not solve a classical SDE, 
while its square does (see Revuz and Yor [UU] Exercise 1.26 p 451]). We thus reformulate the 
equation in an adequate sense and in such a way it has a unique solution (in law). We also prove 
that this solution is stuck at 0 when x > Stt, while it reaches 0 but escapes instantaneously when 
X e (0,87r). 

1.6. References. Approximating a large particle system by a partial differential equation (for 
deriving the PDE) or a partial differential equation by a large particle system (to compute numeri¬ 
cally the solution of the PDE) is now a classical topic, called propagation of chaos. This notion was 
introduced by Kac m as a step to the rigorous justification of the Boltzmann equation. When 
the interaction is regular, the situation is now well-understood, some important contributions are 
due to McKean [24] , Sznitman [35] , Meleard [25] , Mischler and Mouhot [26] , etc. The main idea is 
that one can generally prove true quantified convergence when the interaction is Lipschitz contin¬ 
uous and tightness/consistency (and true unquantified convergence if the PDE is known to have 
a unique solution) when the coefficients are only continuous. Of course, each PDE is specihc and 
these are only formal rules. 

The case of singular interactions is much more complicated. In dimension one, let us mention 
the works of Bossy-Talay [U] and Jourdain m which concern the viscous Burgers equation and 
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more general scalar conservation laws (where particles interact through the Heaviside function) 
and of Cepa-Lepingle [5] on the very singular Dyson model. 

A model closely related to the one studied in the present paper is the 2c?-vortex model, that 
approximates the vorticity formulaltion of the 2d-incompressible Navier-Stokes equation. The PDE 
is the same as m and the particle system is the same as o, replacing everywhere the kernel K, see 
(ED, by the Biot and Savart kernel x'^/{2 'k\x\^). This kernel is as singular as K, but the interaction 
is of course not attractive, so that the situation is simpler. In particular, there is no blow-up for 
the PDE and Osada m has shown that the particle system is well-posed and that particles do 
never collide. Osada [551 has also proved the (true but unquantified) convergence of the particle 
system to the solution of the PDE when y is sufficiently small (in our notation), and this limitation 
has been recently removed in [8]. The method developed in |8] relies on a control of the Eisher 
information of the law of the particle system provided by the dissipation of its entropy. It has 
been applied to a subcritical Keller-Segel equation by Godinho-Quininao [5], where K is replaced 
by —x/(27r|x|^+“) with some a G (0,1) and to the Landau equation for moderately soft potentials 
in [7]. Let us finally mention the propagation of chaos results for some particle systems with 
deterministic dynamics by Marchioro-Pulvirenti [23] (for the 2(i-Euler equation) by Hauray-Jabin 
[m (for some singular Vlasov equations) and by Jourdain-Reygner |18j (for diagonal hyperbolic 
systems). 

In the above mentioned works, some true convergence is derived. Here, we obtain only a 
tightness/consistency result, but the singularity is really strong and attractive. Concerning the 
Keller-Segel equation, we are not aware of papers dealing with the convergence of the true particle 
system without any cutoff. Stevens [3S] studies a physically more convincing particle system with 
two kinds of particles (for bacteria and chemo-attractant particles). She proves the convergence of 
this particle system when the kernel K is regularized. In [10], Haskovec and Schmeiser also prove 
some results for a regularized kernel of the form Kffx) = —x/[|a;|(|x| -I- e)]. Einally, Godinho- 
Quininao [9] study the case where K is replaced by —a;/(27r|a;|^''‘“) for some a G (0,1). 


1.7. Plan of the paper. In the next section, we prove (0 for a regularized particle system. This 
is the main tool for the proofs of Theorem [5] (existence for the particle system when x S (0, 27r)) 
and Theorem [6|(tightness/consistency as V —>• oo when x G (0, 27r)) given in Section [3l as well as 
for checking Theorem [3 (local or global existence for the particle system in the general case) in 
Section m We establish Proposition |4| (positive probability of collisions) in Section [Sj Section |6| is 
devoted to a detailed study of the case N = 2. Finally, we quickly and formally discuss in Section 
[3 how to build an relevant V-particle system when x ^ 8tt{N — 2)/ {N — 1) and we explain why it 
seems to be a difficult problem. 


2. A REGULARIZED PARTICLE SYSTEM 

Let fo G Pi(R^), X > 0 and iV > 2 be fixed. We consider a family Xq, i = l,...,iV of 
/o-distributed random variables and a family {Bl)t>o, i = l,...,iV of 2-dimensional Brownian 
motions, all these random objects being independent. For e G (0,1), we define the regularized 
version of K as 


—x 

2tt{\xY + £^)' 


( 8 ) 


Ke{x) 
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This kernel is globally Lipschitz continuous, so that the particle system 


(9) 


= X^, + V2Bl + j^^ f K,{Xl’^’- - Xi’^’^ds 

i=i 


is strongly and uniquely well-defined. These particles are furthermore clearly exchangeable. The 
following estimates are crucial for our study. 

Proposition 9. For fo G N > 2 and e G (0,1), we consider the unique solution 

iXl’^’'^)t>o,i=i,...,N to ([9]). 

(i) For all t > 0, E[(l + < {fo, ^/T+W) + 2<. 

(ii) For all a G (0,1), all T > 0, all r] G (0, e], 




E 






< 


{2V2{fo,^l + \x\^)+4V2Ty 


jN - 2)x 

ttN 


E 






Proof. We start with point (i). Using the Ito formula (with (j){x) = (1 -I- whence V4>{x) = 

(1 -I- and A(j){x) = (1 -|- |xp)“^/^(2 -|- |a;p) and taking expectations, we find 


E[(l -k =E[(1 -k 




-ds 




N 


i/i 


0 (1 + |X^'^’d2)3/2 

vl,N,e 

K,{Xy’^ - Xi’^’^)ds 


Uo {i + \xy^’y^y/^ 


By exchangeability and oddness of K^, for j G {2,..., iV} 


E 


Uo {i + \xy’y^y/^ 

pt Yl,N,e 

Li • 


= ^E 


• K^{Xl’^’^ - Xf^Uds 

XhN,e 


2 yjo Ml+|X’'^’"P)l/2 (1 + |Xi'^.d2)l/2 


• K^{Xl'^'^ - Xy^Uds 


This last expectation is non-positive since for x,y G the inequality Ixl"^ -I- |y|^ > 2 |xp| 2 /p 
implies (|a:p(l -|- + |2/P(1 + > (|a:||y|((l -I- -I- (1 -I- |a:p)^/^))^, whence 

{x{l + \y\'^Y/'^ — y{l + \x\'^Y/'^)-{x — y) > 0 and thus (a;(l -I- — ?/(l -I- • (cc —y) > 0. 

Hence 


E[(l -k = E[(l -k \X^yY/‘^] + E 


r‘ 2 + \xl’^’y^ 


-ds 


(1 + |X'’'^’"P)3/2 

as desired. To prove point (ii), we fix a G (0,1) and start from 

^l,N,e _ ^2,N,e ^ _ ^2 ^ - B^) + xR^ + XS, 


< E[(l + + 2T 


12 
t 5 


where Ej= 3 /o [^e(X^’^’" - X’^’^) - A"e(X’^’'^ - X’^’'^)]ds and where = 

2N~^ JL Kg{Xy^’^ — Xy^Uds. We next fix 77 G (0,e], introduce X(^) = (kP + and 
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use the Ito formula to write 


=E[ct>r,{X^ - X^)] +e\ I 2Acj,r,{Xy’^ - xy’^)ds 

pT 


+ xE 


- X^’n ■ (dRl^ + dSl^) 


‘-^0 


Since 77 G (0,1), we have (j)Ti{x — y) < [\/2((l + + (1 + \y?'Y^'^))]°‘^ whence £[</>,, 

< (2-\/2(/o, y'l + |a;p) +4-\/2T)“ by (i). Since furthermore E[^^(Xq — Xq)] > 0, 


(10) E / 2A^^(Xi’^’'^ - Xf>^>^)ds 


<(2y2(/o,x/rT^)+4y2T)“-xE / V- X^^^■ {dR^ + dS]^) 


'-Jo 


Using exchangeability and recalling the definition of Rl 
pT 


12 


-E 


V<(.,(Xi’^’^ - . dRf 


^Jo 


N 




But V(j)rj{x) = a(|a:p ^x, whence |V(^^(x)| < Q;(|a:p+77^)“/^ Furthermore, |Xg(a;)| < 
(|a;p + e^)“^/^/(27r). Hence 


(11) -E / V</>,(Xi'^'^-X2X.e)-di?: 


12 




(X — 2)a„^ 
ttN 


(IXlX.e _ ;x2.A,e|2 ^ ^2)a/2-l/2(|;^l.A.e _ ^ e2)-l/2^^' 


'-^0 


Recalling the definition of and using that |X£(a:)| < (jeep + 77^) ^/^/(27r), 


(12 ) -e[ / xuxl’^’^-xy’^)-dsl 


(X „ 

< 

ttN 


(|X]X.e _ ^ 2 ^a/ 2 -l^^ 


Finally, we observe that A0^(x) = a(|a;p + 77^)“/^ ^(a|a;p + 277^) > Q;^(|a;p + 77^)“/^ Inserting 
this into (ITOl) and using (fTT|) and (IT^ . we find 


(13) 2a2E[ / (|XlX.e_^2.A.e|2_^^2^a/2-l^g 


'-^0 


<(2V2(/o, Vl + |*|2) + 4y2r)“ + ^E / (|X]X,^_x 2 X.= | 2 +^ 2 )a/ 2 -l^g 


ttN 


^JO 


{N - 2)a x 
ttN 


E 


(IXiX.e _ x2.A.e|2 ^ ^2)(a-l)/2(|^l.iV.e _ ^3.iV.e|2 ^ £2)-l/2^^ 


'-v'O 


The conclusion immediately follows. 


□ 
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3. Tightness and consistency in the (very) subcritical case 

The aim of this section is to prove Theorems [S] and IH] In the whole section, /q G Pi(R^) is fixed. 
First, we deduce from Proposition |9] an estimate saying that in some sense, particles do not meet 
too much, uniformly in iV > 2 and e G (0,1) when y < 27r. 


Corollary 10. For each N > 2, each x S {0,27rN/{N — 1)) and each e G (0,1), consider the 


unique solution (Xj’^’‘^)t>o,i=i,..,,Ar to For all T > 0 and all a G (x(fV — l)/(27rA^), 1), 


E 


1^. 


l,N.e 


-X 


2,N,e\Oi — 2 


ds 


'-Jo 


< 


i2V2{fo,^l + \x\‘^)+AV2TY 

a{2a -{N - l)x/(7rX)) 


Proof. We thus fix a G {x{^ ~ l)/(2'TiV), 1). By Holder’s inequality and exchangeability, we have, 
for any ry G (0,e], 


Applying Proposition [^(ii), we thus find 

(2a - < 

It suffices to let ry \ 0 to complete the proof. 


(2V2(/o,^l + |a:|2)+4y2r)“ 
a 


□ 


Such an estimate easily implies tightness. 

Lemma 11. For each N >2, each e G (0,1), consider the unique solution (X(’^’’^)(g[Q 
to (O. 

(i) For N > 2 fixed, if X < 2 ttN/{N — 1), the family e G (0,1)} is tight in 

C'([0,oo),IR2). 

(a) If X < 271, the family {{Xl’^’'^)t>o, N > 2, £ G (0,1)} is tight in (^([O, oo), R^). 


Proof. We first prove (ii) and thus suppose that y < 2tt. Since (^([O, oo),R^) is endowed with the 
topology of the uniform convergence on compact time intervals, it suffices to prove that for all 
T > 0, {(X(^’^’^)tg[o_T], X > 2, £ G (0,1)1 is tight in C([0,T],R^). Let thus T > 0 be fixed and 
recall that Xl’^’’^ = Xq + V^Bl + , where 


tI.V.e 

Jt 


X_ 

N 



)ds. 


Observing that the laws of Xg and (Hj )tg[o,T] do not depend on X > 2 nor on £ > 0, it suffices 
to prove that the family {iJt’^’^)t&[o,T], X > 2, £ G (0,1)}, is tight in (^([O,T],R^). To do so, we 
fix a G (x/(27r), 1), and we use Holder’s inequality to write, for 0 < s < t < T, 




<- 


-2ttN 



I ^du 


<\t - s|(i-“)/(2-«)^ ^ (^ ixi’X. _ xi’^’T-^du) 


l/(2-a) 


7=2 
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where (3 = (1 —a)/(2 —a) > 0 and where := (x/(27riV)) + /o^ ~'^du]. 

Indeed, < 1 + a; because a G (0,1). But we immediately deduce from Corollary [TUI and 

exchangeability that sup^g^Qj) jv >2 so that there is a constant Ct, not depending 

on £ G (0,1) nor on iV > 2 such that for all A > 0, P(Zy’® > ^) < CrjA. Since = 0 

a.s., we conclude that for all A > 0, for all TV > 2, all e G (0,1), P[( j/’^’®)tg[o,T] ^ ^a] < CrjA, 
where /Ca is the set of all functions 7 : [0, T] i-A such that 7(0) = 0 and for all 0 < s < t < T, 
|7(0~7('S)I ^ 2 l|t—s|^. The Ascoli theorem ensures us that Ka is a compact subset of (^([O, T], R^) 
for all A > 0. Since limy 4 _j.oo sup^> 2 .£e(o,i) lP[('^/’^’^)tG[o,T] 4- ^a] = 0, the proof of (ii) is complete. 

The proof of (i) is exactly the same: the only difference is that N is fixed so that we can choose 
aG (x(iV-l)/(27riV),l). □ 


We now prove the existence of the particle system without cutoff in the very subcritical case. 


Proof of Theorem[^ We divide the proof in two steps. Recall that y < 2ttN/{N — 1). 

Step 1. For each e G (0,1), we consider the unique solution (Wj’^’^)tg[o_oo).i=i,....Ar to dH). 
By Lemma [TTl-fil. we know that the family e G (0,1)} is tight in (^([O, 00 ), R^). 

By exchangeability, we of course deduce that ,..., ^’®)t>o, e G (0,1)1 is tight in 

(^([O, 00 ), (R^)^) and consequently that , B^),..., {X^’^’‘^, B^))t>o, £ G (0,1)} is tight 

in C'([0,oo), (R^ X R^)^) (this last assertion only uses that the law of {B^,... ,Bf^)t>o does not 
depend on £). It is thus possible to find a decreasing sequence £fc \ 0 such that the fam¬ 
ily , B^),..., , Bf^))t>o converges in law in C'([0, 00 ), (R^ x R^)^) as fc —>• 00 . 

By the Skorokhod representation theorem, we can realize this convergence almost surely. All 


this shows that we can find, for each fe > 1, a solution (X, 


IW.eic 


...,x 


N,N,ek 


)t>o to (I5|), as¬ 


sociated to some Brownian motions {bI’^’^'‘ ,..., B^’^’‘^'°)t>o, in such a way that the sequence 
((X^^’^^ •. ■, a.s. goes to some ((X^^, B }),..., (xf'^, Sf ))t>o 


in (^([O, 00 ), I 


N 


) as fc —>■ 00 . Let us observe at once that the family {(X(’^)t>o, i = 


1,... ,X} is exchangeable and that, by Corollary [10] and the Fatou Lemma, for all T > 0 and 
aG (x(X-l)/(27rX),l), 


(14) 


max ‘ 


|e[ [ |X}’^-X2’^|“-2dsl,supE[ [ 
t LJo J fe LJo -I 


^ {2V2{fo,^l + \x\^)+4V2Tr 

~ a{2a -{N - l)x/( ttX)) 

Step 2. We introduce J) = cr((X*’^, Of course, iXl’^)i=i^,,,^N,t>o is {d^t)t>o- 

adapted. The family (Xq is of course i.i.d. and /o-distributed (because this is the 

case of (Xo’^’®'“)i=i,...,Ar for all k > 1) and (B®)j=i 7 v,sG[o,t] is obviously a 2X-dimensional 
Brownian motion (because this is the case of (i?s^’'^*')i=i,...,Ar,s>o for all k > 1). We now 
show that (B*)i^i_ _jv,sG[o,t] is a 2X-dimensional (Xt)t>o-Brownian motion. Let thus t > 0, 
4> : (^([O, 00 ), (R^)^) and if : C([0, t], (R^ x R^)^) be continuous and bounded. We have to check 
that 


E[^((X:>^, i?:)i=i.....^..e[o.t])</'((SL. - A).=i....X..>o))] 

= E[v;((x:’^, i?:'^),=i.....^.,g[o,t])]i?[<^((sL. - A)i=i....x.s>o))]. 
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This immediately follows from the fact that for all /c > 1, 




i,N,ek 

t+s 


- B 


i,N,ek 


)i=l,...,A,s>o))] 




i,N,Sk 

t-\-s 


- B 


i,N,ek 


)i=l....,A,s>o))], 


which holds true because (-^t’^’^‘’)i=i,...,A,t>o is a strong solution to (0 and is thus adapted to 
the filtration = (7((X^, Bl’^•^'‘)i=l^,,,^N,s>o)■ 

Step 3. It only remains to check that for each i € each f > 0, Xl’^ = X^^ + V^Bl + 


ri,N,ek _ Yi,N,Ek 


{x/^)J2j=iIo^(^s’^ — X^X'jdg^ yve of course start from the identity = X^ 

\/ 2 Bl’^’^k _|_ {x/N)J 2 ^=i Jo Xek{XlX,sk _ and pass to the limit as fc —>■ oo, e.g. in 

probability. The only difficulty is to prove that Jl^ (t) tends to where 

J^(t)= f - Xl’^’->^)ds and = f K{Xl’^ - Xl’^)ds. 

Jo Jo 

We introduce, for r] G (0,1), 


I and {t) = j K^{Xl’^ - XiX)ds. 


For a G (0,1) and k sufficiently large so that Sk < d, we have 


(15) 


\K^{x) - iFeJx)| + \Kr,{x) - K{x)\< 


V 


7r|a;|(|a;p + 


< 


a|2;|“ 2 


We thus deduce from (fHl) that, for a G (x(X — l)/(27rA^), 1), there exists Ca,t < +oo such that 


E[|J*^(t) - j;^(t)|] +lim^supE[|4^(t) - < Ca,tV 


l-a 


Next, since is continuous and bounded and since (Wg’^’®'“)s>o goes a.s. to (Wg’^)s>o, it holds 
that Jlfit) a.s. and in for each ij > 0. Writing 

- Ji^m < E[\r^{t) - j;^(t)|] + E[|j;^(t) - rl{t)\] + E[\rdjt) - Ji^m, 

we conclude that limsupfc_^oo (0 ~ ii)\] ^ Ca,td^~°‘- Since rj G (0,1) can be chosen arbi¬ 
trarily small, we deduce that indeed, {t) tends to in as fc —>■ oo. □ 


Following some ideas of [H Proposition 6.1], we now give the 

Proof of Theorem\^ For each N > 2, we consider the particle system built 

in Theorem O and we set which a.s. belongs to 7^(C'([0, oo),M^)). 

For each t > 0, we also set = N~'^ S^i.N, which a.s. belongs to 7^(IR^). 

Step 1. For each N > 2, {Xl'^)f^[o,oc),i=i....,N has been obtained as a limit point (in law), of 
(X(’^’®)tg[o_oo),i=i....,A as e —>■ 0. By Lemma [Til fiib the family {{Xl’^)t>o, iV > 2} is thus tight 
in C([0, oo), R^). As is well-known, see Sznitman [35l Proposition 2.2], this implies that the family 
{fi^, A^ > 2} is tight in 7^(C'([0,oo),R^)) (because for each N > 2, the system is exchangeable). 
This proves point (i). 

Step 2. We now consider a (not relabelled for notational simplicity) subsequence of going 
in law to some fi and show that fi a.s. belongs to 5 := {£((Ai)t>o) : (At)t>o solution to the 
nonlinear SDE (I3|) with initial law /o}, recall Definition [2| This will prove point (ii). 
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Step 2.1. Consider the identity map 7 = ( 7 t)t>o : 00 ),i-A 00 ),M^). Using the 

classical theory of martingale problems, we realize that Q € 7^(C’([0, 00 ), K^))) belongs to S as 
soon as, setting Qt = Qo 7 ^“^ £ P(R^) for each t > 0, 

(a) Qo = /o; 

(t*) So /r 2 /r 2 | 2 ^ “ y\~^Qs{dy)Qs{dx)ds < 00 for all T > 0; 

(c) for all 0 < ti < • • • < tfc < s < t, all ipi,... ,'pk & all ip G C'^(K^), 

J^iQ) JJ Q{dz)Q{dz)ipi{zt^)...ipk{zt^) 

ip{zt) - ip{Zs) -X dC{Zu - Zu) ■ Vip{zu)du - 


A(p{zu)du 


= 0 . 


Indeed, let {Xt)t>o be Q-distributed, so that C{Xt) = Qt for all t > 0. Then (a) says that Xq is 
/o-distributed, (b) is nothing but the requirement (a) of Definition [2l and (c) tells us that for all 
</pGC2(K2), 


K{X, - Zs) ■Xp{Xs)Qidz)ds - [ Aif{Xs)ds 

JO 

is a martingale in the filtration {J-t)t>o generated by {Xt)t>o- This classically implies the existence 
of a 2-dimensional (J't)t>o-Brownian motion {Bt)t>o such that Xt = Xo + V2Bt + X fo f — 

Zs)Qs{dz)ds for all t > 0. It finally suffices to observe that for all a; G and all s > 0, f K(x — 
Zs)Q(dz) = {K-kQs){x). 

We now prove that y a.s. satisfies these three points. For each t > 0, we set p.t = p.o 

Step 2.2. Since is the empirical measure of N i.i.d. /o-distributed random variables and 
since po is the limit (in law) of Pq , we obviously have that po = fo a.s., i.e. p a.s. satisfies (a). 

Step 2.3. Using Corollary [TUI and exchangeability, we see that for any a G (x/(27r),l), any 
T > 0, there is a finite constant Ca,T such that for all m > 0, all > 2, 


(p(Xt) - (p(Xs) -X J j 


E 



(TOA|a;-y|“ ^)p‘f {dy)p^ {dx)ds 


<■ 


mT 


1 


N m 


i¥=3 


- XI 


,N\Ci-2 


ds 




-lC 


as 


Since goes in law to the LHS converges to E /i^ 2 /i[^ 2 (^ A ~ y\^~^)fJ's{dy)fis{dx)ds 
N ^ 00 . Letting m increase to infinity and using the monotone convergence theorem, we find that 

pT 


E 




'^ps{dy)ps{dx)ds 


< Ca,T- 


Since a < 1, this of course implies that p a.s. satisfies (b). 

Step 2 . 4 . From now on, we consider some fixed X : 00 ), R^)) >->• R as in point (c) and 

we check that X{p) = 0 a.s. 


Step 2 . 4 . 1 . Here we prove that for all N >2, 


(16) 


E 


iX{p^)f 


< 


N ■ 
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To this end, we recall that G C^(R^) is fixed and we apply the ltd formula to (|1]); 

Of (<) ^ E . K(x:’^ - x^’^)ds - f x^{xix)ds 

j=l ''0 Jo 


=<f{X^) + V2 fxifiiXlX) . dBl. 

Jo 


By definition of X (recall that K{0) =0 by convention), 


N 


=N ^ - of (5)] 

f2. 


2=1 

N 


[ x^{xiX).dBi. 

i=l 


Then (fTBl) follows from some classical stochastic calculus argument, using that 0 < ti < • • • < 
tk < s < t, that (fii,... ,(pk,X(p are bounded and that the Brownian motions B^,..., B^ are 
independent. 


Step 2.4-2. Next we introduce, for 77 G (0,1), J), defined as T with K replaced by the smooth 
and bounded kernel recall Then one easily checks that Q 1 —>■ (Q) is continuous and 
bounded from ^((^([0, 00 ), M^)) to R. Since goes in law to p, we deduce that for any p G (0,1), 

E[|J-f/r)|]=limE[|J-f/r^)|]. 

Step 2-4-3. We now prove that for all > 2, all 77 G (0,1), all a G (x/(27r), 1), 

E[\X{p) - X^{p)\] + sup E[| - X^{p^)\] < 77 !-“. 


Using that all the functions (including the derivatives) involved in X are bounded and that we 
have |N'T;(a;) — K{x)\ < /tt by (ITSl) . we get the existence of a finite constant Cjr 

such that 


\X{Q)-X^m <CrP 


1 -a 


I 


\a-2 


l{z^^z^}du Q{dz)Q(dz) 


--Cjrp^ “ / [ / lx-7/1“ '^l{^^y}Qu{dy)Qu{dx)du. 

Jo JR2JB2 


The conclusion then follows from Step 2.3. combined with the estimate 


E 


|x - 7/1“ l{x^y}Ps {dy)p^ {dx)ds 


to Jv?Jv? 

deduced from Corollary [10] and exchangeability. 


< 


1 

w 


Ee[ / \xy^-x^ 


N \ol — 2 


ds 






< Ca,T 


Step 2.4-4- For any 77 G (0,1), we write 

E[|J-(/ 7 )|] <E[|J-(/ 7 ) - X,{p)\] + limsnp |E[| / 7 )|] - E[| J-,( 77 ^)|]| 

N 

+ limsupE[|J^^(/ 7 ‘^) — X{p^)\] + limsupE[|J^(/ 7 '^)|]. 

N N 

By Steps 2.4.1 and 2.4.2, the fourth and second terms on the right-hand side are zero. We thus 
deduce from Step 2.4.3 that E[|J^(/ 7 )|] < Cq^jf 17 ^““. Since 77 G (0,1) can be chosen arbitrarily 
small, we conclude that E[|J^(/ 7 )|] = 0, whence X{p) = 0 a.s. as desired. 
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Step 3. It only remains to check point (iii). Consider the (not relabelled) subsequence going 
to p in V{C{[0, oo),IR^)) as in Step 2. This implies that {p^)t>o goes to {pt)t>o in ^”([0, oo), 7^(]R^)). 
By Step 2, p is a.s. the law of a solution to the nonlinear SDE ([3]). As seen in Remark |31 this 
implies that a.s., {pt)t>o is a weak solution to the Keller-Segel equation o. □ 

4. (Local) existence for the particle system in the general case 

The aim of this section is to prove Theorem!?] We thus fix x > 0 and /o G Pi(K^). We introduce 
the domain, for £ > 1 and N > 2, 

Df := {(ail,..., xat) G (M^)'^ : + a;fc| + |a;fe—a;i| > l/£ for alH,j, fe pairwise different}, 

and we consider the Lipschitz continuous function : (R^)'^ i—>■ [0,1] defined by 

$^(a;i,... ,a;Ar) = 0 V (2£ min {\xi - Xj\ + \xj - Xk\ + \xk - Xi\} - 1] Al, 

\ i,j,k distinct / 

which satisfies Inw < As usual, the random variables Ain, i = are i.i.d. 




with common law /o and independent of the i.i.d. 2-dimensional Brownian motions (i?))t>o, 
i = 1,..., iV. For e: G (0,1) and £> 1, the particle system 

N 

(17) = X* + V2BI + f E / ((X^^'^’')fe=i„„.^)ds 

is strongly well-posed, since and are bounded and Lipschitz continuous. For a fixed i> 1, 
we can show as in Corollary ITOl that particles do not meet too often. 

Lemma 12. Fix y > 0 and consider, for each N > 2, e € (0,1) and £ > 1, the unique solution 

to ini). 

(i) For all t>0, all £> 0, all e G (0,1), E[(l -b ^1-b |a;|2) + 2t. 

(ii) For all T > 0, all a € (0,1), all £ > 0, there is a constant CT,a,i (depending also on x and 
fo) such that for all e G (0,1), all N > xlif^ctn), 


E 


1^1,A, 


'-ao 


< 1 + 


CT,a,f. ( 


2a 


X \(“-2)/(l-a) 


ttNJ 


Proof. First, (i) can be checked exactly as Proposition |9]-(i), using only that <i>^ is nonnegative 
and does break the exchangeability. We now prove (ii) and thus fix a G (0,1). Proceeding exactly 
as in the proof of Proposition |9]-(ii), see (flSl) . we hnd that for all rj G (0,e], 




AQ: jN,e,l (-^ — ‘^)x<x j^N,E,e 




ttN 




where = (2-\/2(/o, \/l -b \x\'^) + 4-\/2T)“ and where 


=E| [ (IM'"'''' - 

^(1^1,A,M ^ ^2yl/2 


jN^€,i _ TT7 


=e 


'-Jo 


'-Jo 
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Since < 1 , we obviously have < I^aT- We next note that for > 0 , < 

(1 + < (1 + max{it, and that, for x = (xi,..., xn) S 

(R^)^, $^(x) > 0 implies that x G whence |xi — X2I + |xi — X3I + |x2 — X3I > 1/(2^) and thus 
max{|xi — X2I, |xi — X3I} > l/(8£). Consequently, since 77 G ( 0 ,e], 

(|a;i - X2|^ + r7^)^““^^'^^(|xi - xsj^ + (x) 

<[1 +max{?7^ + |xi - X2p,r7^ + |xi - X3p}“^/^] 

X [1 + (772 + |xi - X2n"^ + (772 + |xi - X 3 \^)-^] 1 {^^D^^} 

<(1 + 8^)[1 + (|xi - X2p + 772)-1/2 + (1^, _ 3;3|2 + 


This implies that 




(1 + 8^)E J + + + 


<(l + 8£)r + 2(l + 8£)Ef / (|X,iW.e.^_x2,Ar.M|2_^^2)-i/2^^ 


‘-^0 


<(1 + 8 i)T + 2(1 + 8 £) 7 ’(l-«)/( 2 -«)[/^^M]l/( 2 -a) 


by the Holder inequality. All in all, we have checked that 

( 2 « - 

where Ba,T,i = Aa,Tla + (1 + 8 £)Tx/'^ and Ca,T,t = 2(1 + 8 i)T^^~°‘'^x!■ 
Separating the cases ^ easily conclude that 

TN,e4<.,(r. ,r X 

L.a,T ^ 1 + [Ba,T,i + Ca,T,eJ [^2a - — J 

It finally suffices to let 77 \ 0 to conclude the proof. 


□ 


We now deduce some compactness, still for £ fixed. 

Lemma 13. Fix x > 0 and consider, for each A^ > 2, e G (0,1) and £ > \, the unique solution 
(M’^’^’^)t>o.i=i,...,Ar to dlT]). For all £ > 1, the family N > max{2, x/(27r)}, e G 

(0,1)} is tight in C([0, 00 ), I 


Proof. We fix £ > 1 and T > 0. As in the proof of Lemma [TTl the only difficulty is to prove 
that the family ^ > -^O) e S (0,1) > 0} is tight in (^([O,T],R^), where 

No = [max{2, x/(27r)J + 1 and 

N i-t 

1=1 do 

We consider a G (0,1) such that 2a — x/i'^No) > 0, so that, by Lemma [T^ 


(18) 


sup E 
Af>7Vo.Ee(0.1) 




a-2 


ds 


< 00 . 
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Using that < 1, we check as in the proof of Lemma ITT] that for all 0 < s < t < T, we have 

— s|^, where ^ = (1 — a)/(2 — a) and where 

N 


7N,e/ _ X 
^ 2t:N 


i=2 


^ds 


But (fT51) and exchangeability imply that sup^>jY^_^g(g ] < oo. We conclude exactly as 

in the proof of Lemma [TT] □ 


We now make e tend to 0 in the particle system dig, simultaneously for all £ > 1. 


Lemma 14. Let y > 0, fV > max{2, x/(27r)} and fo € 7^i(]R^) be fixed. There exists, on some 
probability space endowed with some filtration {iFt)t>o, a family (Xg)i=i_...^jv ofi.i.d. f^-distributed 
ipQ-measurable random variables, a 2N-dimensional {Tt)t>o-Brownian motion (-Bj)i=i_,,,_ 7 v_t>o and, 
for each i> 1, an {Tt)t>o-o,dapted solution to 


N .t 


(19) 


=Xl + V2Bl + j-^j 


The family i = 1,... ,N} is furthermore exchangeable. Moreover, for all ^ > 1, 

all t > 0, we have E[(l + < (/o, -^1 + |a;P) + 2t and, for all a € (x/i^irN), 1), 


E 


ixi 


l,N,e _ ^2,N,l\a-2 


] < OO. 


Finally, we have the following compatibility property: for all 1, a.s., (X^= 


(X»X>^ for all t € [0,r^), where 


t^n = inf{i > 0 : {Xf^'%=i,...,N i 


Proof. We thus fix y > 0, X > max{2, y/(27r)} and fo G Pi(R^) and divide the proof in several 
steps. 

Step 1. We know from Lemma [13] that for each £ > 1, the family {(X^e € (0,1)} 
is tight in (^([O, oo), R^). By exchangeability, {(X}’^’'^’^)t>o,i=i,...,A, e G (0,1)} is thus tight in 
(^([O, oo), (R^)^), still for each £>1. Since (^([O, oo), (R^)^) endowed with the topology of uniform 
convergence on compact subsets of [0, oo) is a Polish space, by the Prokhorov theorem, for all p > 0, 
we can find a compact subset of (^([O, oo), (R^)^)) such that sup,,g(Q 2 ) P((X}’^’®’^)t>o_j=i_,,,_Af ^ 
X^) < r]2~^. We now introduce X,, := ni>i is a compact subset of [(^([O, 00 ), (R^)^)]^ 

(endowed with the product topology) by Tychonoff’s theorem. It holds that 

sup P(((X*’'^’"’^)t>o.fcl,...,A)^>l ^ x^) < ^ sup P((X*’'^’"’^)t>o.i=l.....A i x^) < p. 

eG(0,l) £>1 

Consequently, the family {((X}’^’‘"’^)t>o,i=i,...,N)£>i, £ S (0,1)} is tight in [(^([0, 00 ), (R^)^)]". 
Finally, we conclude that the family 

{(((X}’^’®’^)t>0,i=l,,,.,A)£>l, (St)t>0,i=l....,A), £& (0,1)} 

is tight in [C([0, 00 ), (R^)^)]” x (^([0, 00 ), (R^)^). 

Step 2. We now use the Skorokhod representation theorem: we can find a sequence £k \ 0 
and a sequence (((X}’^’®'“’^)t>o,i=i,....A)^>i, {Bl'^)t>o,i=i,...,N) going a.s. in [(^([O, 00 ), (R^)^)]^ x 
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C([0,oo, 


N 


) to some {Bl)t>o,i=i,..,,N) and such that, for each £ > 1, 

i,N,ek,^\ _ /TT^ n _:__ j.: _ 


each k > 1, )t>o,i=i,...,A solves (fT71) with the Brownian motions {Bl' )t>Q,i=i,...,N and 

some i.i.d. /g-distributed initial conditions {Xq (not depending on £ > 1). The 

exchangeability of {{Xl’^’^)t>o/>i, i = 1,..., N} is inherited from that of i = 

1,..., N}. Next, Lemma IT^ and the Fatou Lemma imply that for all t > 0, all £ > 1, 

max{E[(l + )i/2],supE[(l + 

and that, for all a G (x/(27r7V), 1), all T > 0, all £ > 1, 


E 


1 ^. 


l,Nl 


-X 


2,Nl\a-2 


ds 


■ sup 

k>l 




< OO. 


Step 3. We introduce Tt = to which is of 

course adapted for each £ > 1. We clearly have X^^’^ = X^^’^ for all i = l,...,iV and all 
^ 1 (because = Xg for all fc > 1, alH = 1,..., and all £,£' > 1). We thus 

may define X^ := X^^’^ for all * = 1,..., A^, for any value of £. The family (Xg)j=i^ consists of 
i.i.d. /o-distributed random variables (because it is the limit of such objects). Finally, one checks 
as in the proof of Theorem [S]- Step 2 (i?j)t>o,i=i,,,,,Ar is 2A^-dimensional (Jb)t> 0 "Brownian motion. 

Step 4- It is checked exactly as in the proof of Theorem [S]-Step 3 that for each £ > 1, 


(X(’^’^)t>o,i=i,,,,,A solves (fT^ : it suffices to pass to the limit in probability as fc —?► oo in the 
equation satisfied by {Xl 
continuous. 


i,N,ehd 


)i>o,i=i,..., 7 V, using the estimates proved in Step 2 and that is 


Step 5. It only remains to prove the compatibility property. We introduce, for £ > 1 and fc > I, 


e.,k 


:= inf{£ > 0 : (X, 


i,N,ekd\ 


^ D^} and := inf{£ > 0 : {XY^'%=i,...,n i 


i,Nd\ 




Since (-^i’^’^''’^)i>o.i=i,..., 7 V goes a.s. to (Ar(’^’®''’^)t>o,i=i,...,A hr (^([O, oo), (M^)^) and since {DfY 


^, we deduce that < liminffc^oo But for all £' > £ > 1, we have 


is an closed subset of , we aeauce tnat S iimmtfc^oo t 
-ri,N,ek/\, , _ _ /Yi,N,ek,i\, , _ ^ In 


[Xl’ )j=i,...,Ar = {Xl’ )j=i,...,Ar on the time interval [0 ,t^ ] for any fc > I: this follows 
from the pathwise uniqueness for (IT^ and from the fact that = 1 on Df. Using finally 

N = {Xl'^'^ )i=i,...,N on [ 0 ,T^] 


that goes a.s. to {Xl’^Xl’^ )t>o,i=i,...,N in C'([0,oo) ''in>2iAf 


^), we conclude that indeed, {X^ 
Finally, we let £ increase to infinity. 


i,N,e\ 


Proof of Theorem^^ We fix y > 0, A^ > max{2, x/(27r)} and /o G Pi(R^) such that fo{{x}) = 0 for 
all X G We consider the objects built in Lemma [T4l the filtration (Jb)t>o, the 2A^-dimensional 
(-^i)t>o-Brownian motion (St)i=i,...,Ar,t>o, the (J'tjt^o-adapted solution (Xt’^’^)t>o,i=i,...,Ar, for 
each £ > 1, to (El), and associated stopping times r^. Using the compatibility property, we 
deduce that is a.s. increasing (as a function of £) and we define tn = sup |>2 r^. Still using 
the compatibility property, we deduce that for all t G [0, tat), all £ such that rf^ > t, all £' > 
£, (Wj’^’^)i=i_...^ 7 v = )i=i.....A- Hence for t G [0,r7v), we can define as 

for any choice of £ such that > t. Since = I for £ G [0, r^], 

by the definitions of and of r^, we conclude that indeed, (^i’^)t6[o,riv).i=i,....Af solves (|1]) with 
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the Brownian motions and that = inf{t > 0 : (>^(’^)i=i,...,Ar ^ Df}, so that 

tn = sup inf {t > 0 : ^ D^} 

i>i 

as in the statement. The exchangeability and (J^t)t>o-adaptation of the family {(^t’'^)te[o,rjv)! * = 
1,... ,N} is of course inherited from i = 1,... ,N}. We also have a.s., for all 

t e [0,Tjv), all a e (x/(27riV), 1), 

f \xi’^ - = f - xf^^^r-^ds 

Jo Jo 

as soon as £ is large enough so that > t. This last quantity is a.s. finite by Lemma [14] again. 

It remains to decide whether tn is finite or infinite. For I C {1,..., N} with cardinality |/| > 2 
and t e [0, r^), let Xj = \I\-^ E.e/ and R( = 2-i E.e/ 1^*’^ - P- 

First assume that x > 8tt{N — 2)/{N — 1). Consider /^r = {1,... ,N}. A direct computation 
using the ltd formula (see (EU)) in the proof of Lemma [T5| below, the last term obviously vanishes 
when / = /jv) shows that )ie[o.TN) is a squared Bessel process with dimension {N — 1)(2 — 
x/47r) < 2, restricted to [0, rjv). But a squared Bessel process with dimension smaller than 2 a.s. 
reaches zero in finite time, see [32l page 442]. We conclude that on the event {rjv = oo}, 
reaches zero in finite time, which of course implies that tjv < oo. Thus P(tjv = oo) = 0 as desired. 

Assume next that x < 8tt{N — 2)/{N — 1). Observe that for {xi,X 2 ,X 3 ) G (R^)^ and x = 
{Xl +X2+ X3)/8, 

\xi - a;2| + |a:2 - - a^il > (|a:i - X2\^ + 1x2 - Xap + jxa - 

= v^dxi - xp + |x2 - xp + |x3 - xp)^/^. 


Consequently, for .^ > 1, 

P(TAr < oo) = P(r7v < oo,T^ < t^) 

i 


= P(< + 00 , . . min inf (|X* - Xl\ + |A,^' - + \X^ - A*]) < )) 

I,J,k distinct t^[0,Ti^) C/ 


< pfxAT <+ 00 , min inf Rj, < -^') ■ 
V 7 : |/|=3tGfO,r„) 61 ^ / 


|7|=3 tG[0,r_iv) 

This last quantity tends to 0 as £ —>■ oo thanks to the following Lemma, whence P(rAr < oo) =0. □ 


Lemma 15. Let N > 8 and x G (0, 87r(A — 2)/{N — 1)]. Consider built 

in the previous proof. For I C {!,..., A} with cardinality j/j > 2 and t G [0,T7v), let Xf = 

andRi = 2-^J2^^j\Xl’^ - X^]^ . For all I C {1,..., N} with |/| > 3, 

PfxTv < OO, inf Af = o') = 0. 

V tG[0,r;v) / 

Proof of Lemma \TS[ We divide the proof in several steps. 

Step 1. Since the initial conditions (AQ^)i<i< 7 v are independent and /o-distributed with 
/o({x}) = 0 for all x G K^, they are a.s. pairwise distinct and for all I C {!,..., A} with 

|/| > 2, P(i?o > 0) = 1. Also, by definition of tjv, we have a.s. Rj > 0 for all t G [0,TAf) and all 
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|/| > 3. For all |/| > 3 and t G [0, rjv), let 


/3t = 


-d= y (X'-” -X!). dBi. 

■Ml ft, 


This process can easily be extended into a one-dimensional Brownian motion (/3/)t>o- In the 
remaining of the step, we check that for t G [0, tat), 

( 20 ) dRl = + (|/| - 1)(2 - ^)dt + ^ -^t)- - Xi’^)dt. 

iti HI 

We work on [0,rAr). Sarting from (|3|) and setting i?/ = Bl, 

d{Xl’^ -Xl) = V2d{B^ - B% + j^[£K{Xl’^ - X}’^) - |/|-i 






dt. 


where Z/ = J2kei ~ ^i' )■ Using the ltd formula, we thus find 

d\Xl’^ - Xi p =2V2{Xl’^ -Xi)- {dBl - dBi) + 

+ -Xi)- [Y,K{Xi^ - Xi'^) - |/|-i 


-iz/ 


Hi 


dt 


and thus 


dRi =V2J2{XI’^ -Xi)- {dBl - dBi) + 2(|/| - l)dt 
iei 

2^ 

N 




iti 


Hi 


dt 


We now observe that - ^i) ■ {dB\ - dBi) = ’ dBj = ^/W^d^i and 

that J2ieiiXt^ — Xi)Zi = 0, so that 

dRi = 2^,dPi + 2(|/| - l)dt +^Y1 

iel Hi 

To conclude the proof of (|20)) . it suffices to note that Xi ■ K{Xi^ — X^’^) = 0 and that 

XlX .K{Xi^ -Xi’^) = \ Y. ■ K{Xl^^ - Xi’^) = - . 

ijeiJH i,Hi,Hi 

Step 2: A key observation. We see in (I20|) that, up to the third-term in the right-hand side, the 
process R^ evolves like the square of a Bessel process of dimension (|/| — 1)(2 — x|/|/(47r7V)). As 
we will show in a few lines, the condition x G (0, SiriN — 2)/{N — 1)] implies that 



(21) min (n — 1)(2 — xn/(47rA^)) > 2. 

n—3,...,N 

Since by [32l page 442] a squared Bessel process of dimension S > 2 a.s. never reaches zero, we 
expect that indeed, for any |/| > 3, R^ a.s. never reaches zero. 


To check (EB), observe that cj){x) = {x — 1)(2 — xa;/(47rA^)) is concave, so that we only have 
to verify that (j){^) ^ 2 and 4>{N) > 2. First, 4>{N) > 2 is equivalent to our condition that 
X < 8'!t{N — 2)/{N — 1). Next, (()(3) > 2 is equivalent to x < 47riV/3. Finally, it is not hard to 
verify that, N > 3 being an integer, we always have 8tt{N — 2)/{N — 1) < 47r7V/3. 
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Step 3. We now prove by backward induction that for all n = 3,..., iV, 

(22) V / C {1,..., N} with |/| = n, < oo, inf = o') = 0. 

V te[o,Tjv) / 

We first observe that (1^^ is clear when n = N. Indeed, |/| = N implies that I = {1,..., N}, 
so that the third term in the right-hand side of (EU)) vanishes and {Ri)te[o,TN) ^ (true) squared 
Bessel process with dimension (TV — 1)(2 — x/(47r)) > 2 restricted to [0,rAr)- Hence inf[o,riv) > 0 
a.s. on the event {tat < oo}. 

We now assume that H22D holds for some n G {4,..., N} and check that it also holds for n — 1. 
We thus consider some fixed I C {1,..., N} with cardinality n — 1. We have to prove that a.s. on 
{tjv < oo}, inf[Q R( > 0. For each j € {1,..., N} \ J, we introduce Ij = lU {}}. 

Step 3.1. We claim that for each j G {!,..., TV} \/, each {xi,... ,xn) G (K^)'^, setting 
x^ = {n- 1)“^ X)ie/ and x^^ = n~^ 

(2n — 3) min \xk — Xj\^ > n ^ \xi — x^^ p — 3(n — 1) ^ \xi — x^\^. 

ie/„ is/ 

We fix fc G / and start from \xk — Xj\'^ = X)ie/1^* ~ ~ i^k 1^* ~ whence, since 

\^i - < 2(n - 2)\xk - + 2J2iei,i^k \^i “ 

(2n - 3)|xfc - Xjp > ^ |xi - Xjp - 2 ^ \xi-Xk\^. 

iGl 

But one easily checks that 2maxfcg/ Y^iapi^k I®* “ - Yi^kai \^i ~ whence 

(2n - 3)mjn|xfc - Xj\^ > ^ |xi - Xj\^ “ X] 1^* “ ^2 ^ ~ ~ 2 ^ ~ 

i,k£l i,kGlj i,kGl 


The claim then follows from the facts that J2ik£i- 1^^ ~ \xi — and that 

Y^.keI\^i -^k? = 2(n - 1) Yiei\^i “ ^^1^- 


Step 3.2. We now fix a > 0 and b = a/3. Step 3.1 implies that when > a and 

2a 


Rj. < b, we have 


(23) 

whence 


min = -, 

keij^i 2n — 3 2n — 3 2n — 3 


\T^/vk,N vj,N^i ^ y/2n — 3 
iiiax |iV(Xt’ -Xi’ )| < ^ . 

keij^i 2'K\/2a 


Hence one may bound the third term in the right-hand side of (EUl) from below: 


(24) 


1 r. ^ 

>a,Rl<b} ]\^ 




> - 


lei j0 

xV2n - 3 
2'kN\/^ 




iei 30 


with c := {N -I- 1 — n)x\/(2ri^-'3y(ri^^T)/(27rTV-ya). Let us now define the stopping time 

CTa = inf {t G [OjTjv) : mini?(^ < a} 

10 
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with convention inf0 = tn and introduce the process (-Rt’°)tg[o,Tiv) defined by = R{ for 
t € [0,(Ta) and, when da < tn, by being the unique solution, for t € [cra,TAr), to 

rt 


pJ _ pi 


+ (|/| - 1)(2 - {t - da). 


The existence of a pathwise unique solution to this equation follows from [3U Theorem 3.5 p 390]. 
We deduce from (001) that this process satisfies, for all t € [0,riv), 

rt 


=R^ + 2 


X 


^d^f + (|/|-l)(2-|§)t 


■ - ^s) ■ K{Xl’^ - Xl’^)ds. 

“'o zG/ j^I 

Step 3.3. Recall that a > 0 and b = a/3 are fixed and that c > 0 has been defined in Step 3.2. 
The existence of a solution (R*’^)t>o such that P(Vt > 0, R^’^ G (Oi^]) = 1 to the SDE reflected 
at the level b 


(25) 


R^’" = Ro^ A 6 + 2 /J + (|/| - 1) (2 - 2^) t - c J* 


-Lt 


(Ls) s>o is an adapted increasing process such that Lq = 0 and (6 — Rl’^)dLs = 0 

will be checked in Step 4 using that |/| > 3. We take this for granted and show that a.s., for all 
tG [0,riv), R^“>R[’^ 

By [32l Lemma 3.3 p 389] with the choice p{u) = juj, the local time at 0 of the continuous 
semimartingale St = Rf' — Rf'^" vanishes. Indeed, it suffices that a.s., /p (p(S's))“^d(S', S')s < 00 , 

which follows from the fact that d{S,S)s = ^[\JR^s^ — \/ Ri’°‘)‘^ds < 4|Rf’^ — Rf’“|ds = Ap{Ss)ds. 

Hence, setting a;+ = max(a;,0), one has, by Tanaka’s formula, for all t G [0,rAr), 

(Ra" - Ri’T =(^0’' - Ro’“)+ + 

Since Rq’^ — Rq’“ = Rg A 5 — Rq < 0, we find 


(Ry^-Rt’“)+< 2 / 1 






- \f^^)dl3i - I' 




ds. 


“'o ie/ 

Since L is an increasing process, the second term on the right-hand side is nonpositive. The third 
term on the right-hand side is also nonpositive, because s < da implies that Rf’“ = R^, so that 
Rg’^ > Rf’“ implies that Rf < b, whence, using (1^ and the definition of da, for all s G [0,tjv) 
such that R{’^ > R{’“, 


-^1 


N 


{s<cra} 


^ J2{Xl’^ - X/) • K{Xl’^ - Xi’^) < c^Rl = c\IrI'^ < cjRi'^ 


iei j0 

We conclude that a.s., for all t G [0,T7v), 


(26) 


(Ra' - R[’“)+ <2^ l{RrA>«L.j(VR?- 
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We next introduce Mt := /q 1{s<tjv y —~ \/ Ri'°‘)d(3l, which is a true martingale 

(because the integrand is clearly bounded by v^), which is a.s. nonnegative for all times by (l26l) 
and which starts from 0: we classically conclude that a.s., Mt vanishes for all t > 0. Coming back 
to (l26l) . we deduce that {rI'^ — Ri’°')~^ < 2Mt = 0 a.s. for all t £ [0, r^r), which ends the step. 

Step 3.4- We now conclude the induction. For any a > 0 and b = a/3, using that {Ri)t£lo,aa) = 
) and the definition of CTq, 

IPItat < oo, inf Rj = 0 

V te[0,r,v) 

= pfrAT < oOjCTa = Tjv, inf = o') +P(^TAr < oo,min inf < a] 

\ tG[0,TN) ' V tG[0,TN) / 

< PI^tat < oo, inf rI’^ = o') +P(^TAr < oo,min inf i?/’ < o'). 

V ie[o.T„) / V tG[o,TN) J 

Since the continuous process (Rt’^)t>o does not reach 0, the first term in the right-hand side is 0. 
We thus can let a tend to 0 to get 

p[t]v < oo, inf = 0 ) < Pfrjv < oo,min inf rI^ = o'). 

V tG[0,r„) J V ie[0.T„) ) 

This last quantity vanishes by our induction assumption. 

Step 4- To conclude the proof, we still have to check the existence of a solution (i?t’^)t>o such 
P(Vt > 0, G = 1 to (1251) . For i > 1/6, according to Skorokhod [33], existence and 

trajectorial uniqueness hold for the reflected (at 6) stochastic differential equation with Lipschitz 
drift and diffusion coefficients 

fA 6 + 2 f* + (|/| - 1 )(2 - ^)t - c/o‘ - Lf 

Wt > 0, < b 

[ {LI)s>o is an adapted increasing process such that Lq = 0 and f^(b — Ri’^'^)dLg = 0. 

Denoting by lye = inf{t > 0 : Rf’^’^ < 1/^}, we deduce from pathwise uniqueness that for £' > £, 
,Tf)te[o...] and (E[’'’', coincide and thus that £ i—>- i^e is a.s. increasing. Setting 

zzoo = we easily deduce the existence of a solution (Rf’^j Rt)tG[o,i^oa) (|25ll satisfying 

sup(g[Q < 6 and Rj’^ > 0 for all t € [0, z^oo). More precisely, Rt’^ = > l/£ for all £ 

and all t £ [0, vi). It thus only remains to prove that Vao = 00 a.s. 

By the Girsanov theorem, under the probability measure Q defined by ^liT(fly(/3r)„g[o t]) = 
exp(c/3//2 — c^t/8) (which is of course a true martingale), the process Wt = jSj — ct/2 is a one¬ 
dimensional Brownian motion. We introduce the equation, satisfied by {rI’^^jL^) on the time- 
interval [ 0 , 1 ^ 00 ), for a squared Bessel process {pt,Rt)t>o of dimension (|/| — 1)(2 — x|/|/(47rfV)) 
driven by W and reflected at the level 6, 

fp, = A 6 + 2/J ^dWs + (|/| - 1)(2 - ^)t- At 
(27) Wt > 0, pt < 6 

I (As)s>o is an adapted increasing process s.t. Aq = 0 and f*{b — ps)dAs = 0. 

To check global existence for this equation, we set Vq = 0 and define, inductively on A: > 0, pt to 
be equal to 


)<P( 


< PI Tat < 00 , CT( 


= ttv, inf 

tG[0,riv 


Ri = 0 


') +P( 


Tat < 00, CTa < Tm 
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• the squared Bessel process 

Rt = 'i-{k=o}Ro A 6 + l{fc>i} 3 ^ y + (|/| — 1) ^2 — ^ ~ 

—k 

on the time interval [ 77 ^, 77 ^+ 1 ] where fjk+i = inf{t > 77 ^ : Rt > 26/3}, 

• the solution to the stochastic differential equation with Lipschitz coefficients 

i?/=I+2 y‘ ~ ’ 

reflected at 6 on the time interval [? 7 fc+i ,?7 ] where 77 = inf{t > 77fc_|_i : R\ < 6/3}. 

Since, under Q, the delays ( 77^+1 — 'rjk)k>i are i.i.d. and positive, Q-a.s., fjk goes to 00 with k by 
the law of large numbers and pt is defined for t € [0,+ 00 ). It is easily checked that the process 
(A)t>o defined by the first equality in also satisfies the last one. 

Reasoning like in the comparison between and performed in Step 3.3, we check that the 
first component of any of two solutions to (1271) is above the other one so that the first components 
coincide. 

We deduce that rI’’’ and pt coincide for t G [0, Uoo). With the definition of Uoo and the continuity 
of p, this implies that {voo ^t} C {3s € [0,t] : ps = 0}. Since {pt)t>o always evolves as a squared 
Bessel process of dimension (|/| — 1)(2 — x|/|/(47rW)) > 2 under the level 6/3, by [32l p 442], 
Q(3s G [0,+oo) : Ps = 0). For each t G [0, 00 ), we deduce that 0 = P(3s G [0,6] : Ps = 0) > 
P(uoo < t) by equivalence of P and Q on cr(i?Q, (/?})sg[o,t]). Letting 6 — 00 , we conclude that 
P(uoo < 00 ) =0. □ 


5. Positive probability of collisions 

The goal of this section is to establish that in the A^-particle system, pairs of particles do collide. 
The main idea is that for e.g. / = {1,2}, up to the third term in the right-hand side of (l20ll . the 
process i?/ resembles a squared Bessel process with dimension (2 — x/(27rA^)) < 2, which a.s. 
reaches 0 by [321 page 442]. 

Proof of Proposition^^ We thus consider any fixed TV > 2, x > 0, /o G P(K^), to > 0 and any 
solution (if it exists) to (|4|). We work by contradiction and assume that a.s., 

Xl’^ ^ for all s G [0, to] and all i ^ j. Then the singularity of K is not visited and the particle 
system (|4|) is classically strongly well-posed on [0,to]. Thus for ff^-a,.e. {x^,... ,x^) G 
there is a unique strong solution jv,te[o,to] to ([4|) such that a.s., XqX = x'’ for all i and 

XlX ^ XiX for all s G [0,to] and all i ^ j. We fix for the rest of the proof an initial condition 
... tX^) G enjoying these properties. All the processes below are defined on the finite 

time interval [0,to]. 

Step 1. By construction, d = mini^j \x’' — x^\ > 0 and we may of course assume that d = \x^ — x'^\. 
We introduce x := [x^ +x^)/2 and note that min 3 <j<Ar \x^ — x\ > 'JZdj^. Fix 1/2 < a < 6 < •\/3/2 
and consider the stopping time r = min{Ti, r 2 , T 3 }, where 

Ti = inf |t G [0,to] : {X^X _ x/’^] > ^ d|, 

r 2 = inf |t G [0,to] : {X^X + - 2^1 > 

T 3 = inf |t G [0, to] : — x\ < bdL^, 
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with the convention that inf 0 = to- We will use that a.s., for all t G [0, r], 

2=1,2, j=3,...,N 

Indeed, consider e.g. the case z = 1 and j = 3, write \Xl'^ — Xf’^\ > ~ ^1 ~ \^t’^ ~ ^1 

and use that -x\>bd and that \Xl’^ - x| < \Xl’^ - X^’^\/2 + {X}'^ + X^'^ - 2 x |/2 < 

(2a + l)d/4 + (2a — l)d/4 = ad. 

Step 2. Consider the exponential martingale defined on [0,to] by 


Mt = exp 


lV2iV^ 


^ pt/\T N 2 

E / E + 1{*>3} E 


N _ ^i,N 


X 

4A^2 


N 


rtAr 


N 


i-1 

2 


E / ib<2} E +1{.>3} E - ^E) 


■dBl 


ds 


i=l ■ 


i=3 


i=i 


This is indeed a true martingale, because K{X^X — Xl’^) is bounded by {2-K{b—a)d)~^ on [0, r] for 
each * = 1, 2 and j = 3,..., TV, see Step 1. Hence P := • P is a probability measure equivalent 

to P. In particular, it also holds that P-a.s., Xl’^ ^ X^X for all s G [0,to] and all i ^ j. The 
Girsanov theorem tells us that, under P, the processes 


IT/ := Bl 


X 


V2N 


ptAr ^ 2 

/ (l{*<2} E + lb>3} E 

•^0 i=3 i=l 


are independent two-dimensional Brownian motions on [0,to]. We next introduce 



(Xb^ - xl^) 

\xix _x^X\ 


■ d 


(W^-WA 
V V2 J 


and 7 t 


Wl + w? 
V 2 


It is easily seen, computing brackets and using Karatzas and Shreve PUI Theorem 4.13 p 179], 
that still under P, /3 is a one-dimensional Brownian motion on [0,to]: 7 , are two- 

dimensional Brownian motions on [0,to]i and all these processes are independent. 

Step 3. We have 


XXN _ ^2,N _ ^2 ^ _ ^ 2 ) + ^ £ K{XlX _ XlX)ds 

+ # E /* {k{XIX - XlX) - K{xy - 

i=3 -^0 

=X^ _ wl) + K{XlX _ xy)ds 

for all t G [0, r]. By the ltd formula, Yt = IX^X _ thus solves, still for t G [0, r], 

T+■ ^y- 
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We also have, for all t £ [0, r] 

N i-t 

=2x + V2{Bl + i?2) + ^ ^ - Xi’^) + K{xy - X^’^)) ds 

=2x + y/2{Wl + W^) 

=2x + 2^t, 


and, for all t S [0, r] and all* = 3,..., iV (recall that K{0) = 0), 




=x* + V2BI + ^ 
=x* + V2WI + ^ 



Xi’^)ds 

-XlX)ds. 


We introduce (yi)te[o.io] ttie unique strong solution, see [32l Theorem 3.5 p 390], to 

y.:=^+2|^dA+(2-^),. 

We clearly have {yt)tG[o,T] = (^t)te[o,T]- We next consider the system 

Xl^ = a;* + ^wi + ^ / E * = 3,..., iV, 

X Jo ^.^3 

which classically has a unique strong solution (X^’ )i= 3 ,,,,,A,tG[o,cr) up to tr = limi^ooinf{t £ 
[ 0 ,to] : inin 3 <i<j<jv — Xp^\ < !/£} (convention : inf 0 = to), which is a.s. positive 

because the initial conditions x^,... ,x^ are pairwise different. Clearly, {X^’ )i= 3 ,...,Ar,tg[o,rAo-) = 
(Xt’^)i= 3 _..._ 7 V.tg[o.rA<T)- We conclude this step mentioning that the processes (Tt)ie[o,to], (7i)te[o.to] 
and (^t’^)i= 3 ,...,Ar,te[o.cr) are independent under P. 


Step 4- For any sq £ (0,to); we claim that 


rii n ^2 c n3 c 



where 


fli = < min Yg = 0 , max W < 
t fo.s ■ ■ ■ 


( 2 a + l^d^ 




03 = { 


= -j (T > Sq, min — xl > bd\. 

se[o,so],i>3 J 


Indeed, on fli, we have max[o,so] < (2a +l)^d^/16, whence, since \Xl’^ — X^X^^ _ 4 ^^ qj^ [ 0 , t ], 
max[o_soAT] < (2a+l)d/2 and thus Ti > sqAt. Since _ 2a ;+274 on [0, r], 

we deduce that on ^ 2 , maxfQ sj,/^,-] + X^X _ 2 a;| < supjo^^^^^] 2 | 7 ,,| < ( 2 a - l)d/ 2 , whence 

T 2 > So A T. On ^ 3 , since cr > so, we have (XjW).^ 3 _ 

and thus minsg[o_soAT],a >3 — x| > bd, so that ro > so A r. As a conclusion, t > sg A t and 
thus r > So on III n O 2 n ^^ 3 - We deduce that Hi O 112 O O 3 C > ^ 0 , minjo^sp] Ts = o| C 

I min[o_so] lA'jW _ = o|, because Yt = jX^X _ for all t g [ 0 ,t]. 
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Step 5. Here we show that we can find sq G (0,to) such that P(ili fl 172 H > 0. As seen at 
the end of Step 3, the events f2i,f22 and Ida are independent (under P). It obviously holds true 
that P(172) > 0 (for any sq > 0) and that P(173) > 0 if sq > 0 is small enough because cr > 0 a.s. 
and by continuity of the sample-paths (at time 0, we have minj >3 — a;| = minj >3 \x^ — a;| > 
\fZdl2 > bd). It thus only remains to verify that P(17i) for all sq S (0, to). Since, by the comparison 
principle stated in [521 Theorem 3.7 p 394], P(ni) is non-decreasing with y, it is enough to check 
that P(17i) > 0 for all sq G (0,to) when x < 47riV, which we now do. 

It holds that T is a squared Bessel process of dimension 6 ■.= 2 — x/(27rAf) started at y = d^/4 
and restricted to the time-interval [0,to]- We set z = (2a -I- l)^ci^/16 and observe that z > y. For 
a; > 0, we also introduce = inf{t G [0, to] : Yt = 2^}- Then Hi = {tq < sq A t^}. 

For a; > 0, we denote by Qx the law of the squared Bessel process of dimension <5 starting from 

X (on the whole time interval [0,oo)), and by qs{x,u) the density of its marginal at time s > 0, 

which is a positive function of u on (0, -boo) according to [521 Corollary 4.1 p441]. For all u ^ v, 

we define Tu as the first passage time at u and Tuv as the first passage time at v after r^. It holds 
that P(r2i) = Qyixo < So A Tz) and what we have to check is that Qy(To < sq A T^) > 0 for all 
So G (0, to). 

We first show that Qx{to < t) > 0 for alH > 0 and all a; > 0. Since 5 < 2, we know from [521 
page 442] that Qx{to < oo) = 1 for all x > 0. With the Markov property, we deduce that 

f+ 00 


1 = ^ Qx(to G (nt/2, {n + 1)7/2]) < Qx^Tq < 7/2) -|- / Q„(ro < 7/2) f ^ qntl 2 {x, u)^ du. 

n>0 '^n>l ^ 

Since, u i—>■ qt/ 2 {x,u) is positive on (0, -boo), this ensures the positivity of 

p+oo 

Qx{to < 7/2) -b l{Q,(ro<i/2)=0} / Qu{to < 7/2)gt/2(x, u)du < Qx{to < 7). 

Jo 

Using the strong Markov property, that 0 < y < z and the monotonicity of 7 >->■ Qy{To < 7), 
Qyi'^z ^ Eo ^ 7) — QyiXzy '^0 — t) — ^{Tzy<t}Qy{,^0 — 7 s)\s—Tzy^Qy — Qyij'zy ^ ^)Qy{^0 — 7)* 

By continuity of the sample-paths, lims-s-o Qyijzy < s) = 0 and we can find si G (0, to) so that for 
all So G (0, si], Qy{Tzy < So) < 1. We conclude that for all so G (0, si], 

Qyiro < So A Tz) = Qyiro < so) - Qyixz < To < So) > (1 - Qy{Tzy < so))Qi/(eo < So) > 0. 

If now So G [si, 7o], we obviously have Qyi^ro < so A r^) > Qy(To < siAtz) > 0. This ends the step. 

Step 6. We deduce from Steps 4 and 5 that P(min[o^tj,] = 0) > 0. But P and P 

being equivalent, this implies that P(min[o_t(,] |A/W _ A/W| = q) > 0, whence a contradiction. □ 

6. Two PARTICLES SYSTEM 

In this section we consider the particle system o with N = 2. Assuming that {X^,Xf)t>o 
solves @ with N = 2, we easily find that St = A/ -b Xf and Dt = A/ — A/ solve two autonomous 
equations, namely St = So + 2Bt and 


( 28 ) 


Dt = Do 


2Wt+x [ K{Ds)ds 

Jo 
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with the two independent 2-diinensional Brownian motions Bt = [B] + Bl)l\f2 and Wt = {B} — 
B^)/\/2. The equation satisfied by {St)t>o being trivial, only the study of (^5)) is interesting. 
During the whole section, the initial condition Dq is only assumed to be a R^-random variable 
indpendent of {Wt)t>o- 

Remark 16. Theorem\^ ensures us existence for (|28l) when y < Itt and Dq is the difference of 
two i.i.d. integrable random vectors. When y > Itt, the equation (1281) has no global (in time) 
solution in the usual sense. More precisely, assume that it has a global solution {Dt)t>Q. Then 
T = inf{t >0 : Dt = 0} is a.s. finite and a.s., \K{Ds)\ds = oo for all h> 0. 

Proof. Let thus y > 47r and assume that there is a global solution {Dt)t>o to (E5)) . By a direct 
application of the ltd formula, this implies that Rt = \Dt\^/A solves Rt = Rq + 2 y^\Rs\df3s + 
(2 - y/(47r))t, where (3t = fg + fg 1 {Ds=o}^A is a 1-dimensional Brownian 

motion (here /3 is any one-dimensional Brownian motion independent of {Dq,W)). According 
to [32l p 442] combined, when y > Stt, with the comparison theorem [32l Theorem 3.7 p 394], 
T = inf{t >0 : = 0} is a.s. finite. By the strong Markov property (for the process R), the 

comparison theorem [32j Theorem 3.7 p 394] and since y > 47r, (RT+t)t>o can be bounded from 
above by a squared 1-dimensional Bessel process starting from 0, process with the same law as 
i\(3tn>o. For h > 0, by the occupation times formula [32l Corollary 1.6 p 224], Jg |/3s| ^ds = 
Ir \a\-^Llda. But > 0 as soon as h > 0 and we know from [32] Corollary 1.8 p 226] that 
a !->■ is a.s. continuous, so that Jg |/3s|“^ds = oo for all h > 0 a.s. Thus 47r |Ar(Zls)|c?s = 
jT+/i all h > 0 a.s. □ 


Hence (1231) has no global solution for y > Itt, while we expect that in some sense, the dynamics 
it represents is meaningful at least for all y £ (0, Stt). We thus would like to refomulate it, in such 
a way that it is possible to build global solutions. More precisely, we would like to identify, for any 
value of y > 0, the limit, as e > 0, of the smoothed equation 

(29) Dt = Do + 2Wt+x f KeiDDds, 

Jo 

where ATg was defined in (jS]). The regularized drift coefficient being Lipschitz, existence and 
trajectorial uniqueness hold for this SDE. We introduce the equation formally satisfied by Zt = 
\Dt\^Dt for {Dt)t>o solution to (125)) : 


(30) 


Zt = Zo+ / <7{Z,)dWs + / b{Zs)ds, 


where o-(z) = 2|2:| ^/^(| 2 :p /2 -I- 22 : 2 :*) and &(z) = (16 — 3y/(27r))|2:| Here and below, /2 is the 

identity matrix and z* is the transpose of z. Here is the main result of this section. 


Theorem 17. Set Zq = \Dq\^Dq. 

(i) 4/y £ (0,87r), (I30l) has a unique (in law) solution {Zt)t>o such a.s., l.{Zt=o}dt = 0. 
Moreover, if x ^ (07 47r), ((28l) has a unique (in law) solution. 

(a) df X ^ Stt, (1301) has a pathwise unique solution frozen when it reaches 0 (and it a.s. reaches 

0;. 

(Hi) In any case, the solution (Ilf)t>o to (1221) goes in law, as £ —>■ 0, to {Dt)t>Q defined by 
Dt = \Zt\~^/^Ztl{Zt^o} and, when y £ (0,47r), this process {Dt)t>o solves (1251) . 
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In point (i), uniqueness in law cannot hold true without restriction for (1301) : the time passed at 
0 by the solution that we consider is Lebesgue-nul, while it is easy to build a solution by freezing 
the process when it reaches 0. 

The rest of the section is devoted to the proof of this theorem. The following lemma is more or 
less standard. 

Lemma 18. Let y > 0 &e fixed. For each e € (0,1), we consider the unique solution (D^)t>o to 
(1^ and we put Zf = |£)f pDf. 

(i) The family {{Zf)t>o,e £ (0,1)} is tight m C^O, oo), R^). 

(a) Any limit point {Zt)t>o is a weak solution to dani) and, setting Rt = \Zt\'^/^/A, it holds that 
(o-) if X € (0, Stt), then {Rt)t>o is a {2 — x/dimensional squared Bessel process; 

(b) if X ^ 8'^) then {Rt)t>o is a (2 — x/dimensional squared Bessel process frozen when it 
reaches 0. 

Proof. We divide the proof in several steps. 

Step 1. Direct applications of the Ito formula show that 

a{Zl)dWs + [ b,{Zl)ds, 

Jo 

where bg{z) = — (3%/(27r))(|zp/^ + £'^)~^z and that i?} := |Df p/4 solves 

«? = «„ + 2 /■ + /' (2 - ;;p^) 

where \{Di=/^o}\Dl\~^Dl‘dWs- Since sup,,>Q(xr)/[27r-y/r(£:^+ 4r)] = x/(87re), the Girsanov 

theorem ensures us that for all T £ (0, +oo), the law of (Rt)te[o,T] is equivalent to the law of the 
restriction to the time interval [0, T] of a 2-dimensional squared Bessel process starting from Rq. By 
[311 p 442], we deduce that a.s., for all t > 0, i?/ > 0. Asa consequence {l3f)t>o is a one-dimensional 
Brownian motion. 

Step 2. By trajectorial uniqueness for (1^ . for M > 0, on the event {|ilo| < Mj, the so¬ 
lution starting from Dq coincides with the one starting from Do1{|£)(,|<m}- Therefore, by both 
implications in the Prokhorov theorem, to check that the family {{Zf)t>o,£ £ (0,1)} is tight in 
(^([O, oo), R^) it is enough to do so when B>o is bounded. The tightness property then easily follows 
from the Kolmogorov criterion, using that sup^g^Q ^^) l^e(^)l and \<j(z)\ both have at most affine 
growth: one classically verifies successively that for all p > 2 and all T > 0 there is Ct,p such that 
for all e € (0,1), sup[o,T] E[|^/P] < Ct,p and E[|Zf - Zf p|] < Crjt - for all 0 < s < t < T. 

Step 3. Using martingale problems, that b and a are continuous and that b^ converges (uni¬ 
formly) to 6, it is checked without difficulty that any limit point {Zt)t>o (as e —>• 0) of the family 
{{Zf)t>o,£ > 0} is indead a (weak) solution to (1301) . 

Step 4- Here we assume that x ^ (0,87r) and we prove that (i?/)t>o goes in law to the squared 
(2 — x/(47r))-dimensional Bessel process. We consider the (2 — x/(47r))-dimensional Bessel process 
{Rt)t>o associated to (/3/)i>o, that is Rt = Ro + ‘2 fg •\/^d/3f-|-(2 —x/(47r))t (its law does of course 
not depend on e) and we prove that lim£_,.o E[sup[o t] l^t ~ ^t|] = 0 for all T > 0, which clearly 
suffices. 

Since /p 3/2 x~^dx = \og{\/e)/2, one may construct a family of nondecreasing convex func¬ 
tions pe : [0,oo) !->■ [0,oo), indexed by e G (0,1/2) such that (pe{x) = 0 for x < '■p'{x) = 1 for 
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x>£ and < C'l{g3/2<2,<£}/[a;log(l/e)] for some constant C € (l,+oo) not depending on e. 

Such functions are called Yamada functions in the literature. We then observe that Rl > Rt for 
alH > 0 by the comparison theorem stated in [3U Theorem 3.7 p 394]. Computing Rl — Rt and 
applying the ltd formula, we obtain that 

^.{Rl - Rt) =2 [ ^liRl - Rs)i^/m - ^/Rs)dPs + ^ 

Jo 471 Jq + ARl 

^0 


We next remark that 


< ^liRl - Rs) 


£2 + 4(i?| - i?,) 




£2 + 4(i?| - i?,) 


< 


V~e 


and that 

- Rs){Vm- < ^e{Rt - Rs){Rl - Rs) < 
whence (the constant C may now change from line to line) 

( 31 ) ^,{Rl-Rt)<2 + 

Taking expectations, we conclude that E[ipe{Rl — Rt)] < Ct/log(l/£). But since (fie{x) < x < 
(Pe{x) + £, we deduce that E[Rl — Rt] < £ + Ct/log(l/£). Coming back to (l3T]) . using the Doob 
inequality and that 0 < < 1 and < Rl — Rs, we conclude that E[sup[Q <fe{Rt — 

.Rt)] ^ C'R/log(l/£) + C^eT + CT^/log(l/£))^/^ and, finally, that E[sup[Q(i?f — i?t)] < £ + 
Cr/log(l/£) + C{eT + C'T2/log(l/£))^/^, from which the conclusion follows. 


Step 5. Finally, we assume that y > Stt and we prove that (i?f )(>o goes in law to the (2—y/(47r))- 
dimensional squared Bessel process frozen when it reaches 0. We consider the frozen (2 — y/(47r))- 
dimensional squared Bessel process associated to (/3( )t>o, that is = i?o + 2 y/Rldfil + (2 — 
y/(47r))t for all t G [0,t], with r = inf{t >0 : = 0} and Rt = 0 for all t>T. We will check 

that for all a > 0, all T > 0, lime_>o P(sup[o,T] |Rt — Rt | > cr) = 0 and this will complete the proof. 
We introduce Tk = inf{t >0 : Rt < 1/fc} and observe that t = supj.>i rj,. 

Step 5.1. For any a > 0, t > 0 and fc > 1, lime_).o P(sup[g ]Rl — i?s| > a) = 0. Indeed, 

using that Rf > Rt for all t > 0 by the comparison theorem [3U Theorem 3.7 p 394], that 

Rl-Rt =2 /Q(i/^-v^)d/3f+ (y/47r) [£2/(£2+4i7J)](is for allt G [0, r^], and that < 

k^/'^]x — y]/2 for all x,y > 1/fc, it is easily checked, by the Doob inequality, that 


E 

sup {Rl - Rs)^ 

< Ck 

[ E 

sup {Rl - Ru)'^ 


[O.tATfe) 


/o 

'-[O.sArfc) 


ds 




whence E[sup[Q t/^T^,) (R| — Rs)^] < Ce'^k'^t^ exp{Ckt) by the Gronwall lemma. 
Step 5.2. We write, for a > 0 and fc > 1 fixed, 

p( sup(i?t - Rt) > a) <]Pf sup (i?t - Rt) > a) + Pf^fc < T, R% > 

^ [O.T] 1 ^ [O.TATfc] 1 ^ 

+ Pfr/c < T, R% < 2/fc, sup Rl> a\. 

^ [ri-.Tl ^ 
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For the last term, we used that — Rt) > a implies that supj^^ Rl > a because 

0 < Rt < Rf- By Step 5.1, the two first terms tend to 0 as £ —>■ 0 (recall that Rr^, = 1/fc), whence 

limsupP^ sup(i?j — Rt) > a') < limsupP^Tfe < T,R^ < 2/k, sup R^ > a). 

p_i.n V jQ jnj / V 1 -^. Ti / 


[rk.T] 


Using the strong Markov property for the process R^ as well as its monotony with respect to its 
initial condition (by the comparison theorem), we deduce that 

limsupP^ sup(i?j — Rt) > C() < limsupP^ sup 

r—^0 V [0,T] ' p—vD V rn T"1 


> a 


• [0,T] 


)■ 


where 


R2jk,e ^2lk + 2 


RV^'^djSl 


ly 


xR. 


2fk,e 


7r(£2 + 4i?r 




We introduce, for r G (0,1) and £ G (0,1/2), the solution (S'/’^)t>o to 5/’® = r + 2 ^J\Sl'^\dj3l + 
2e^ f*(e^ + 4|S'/’'^|)“^ds. Such a solution is pathwise unique by [SH Theorem 3.5 p 390] and 
nonnegative by the comparison theorem [321 Theorem 3.7 p 394]. Again by the comparison theorem, 
and since y > Stt, we find that a.s., for all t > 0. Hence 


limsupP(" sup(i?( — Rt) > a) < limsupP(" sup > a') < limsup 

e^O rn Tl ' p^n rn T”! ' p—^f) 


E[sup[o^T] S\ 


2/fc,£:i 


■ [0,T] 


E —^0 


- [O.T] 


£—^0 


We will verify in the next step that (if r G (0,1]) 
(32) E 


sup < C(1 + T)(t + 1/ log(l/£))l/^ 

[O.T] J 

SO that limsupg^QP(sup[Q y](i?/ — Rt) > a) < C{1 + T)k~^/‘^/a. Letting k tend to infinity, we 
conclude that, as desired, limsup£_jQP(sup[Q 7 ^](i?/ — Rt) > a) = 0. 

Step 5.3. To show (15^ . we consider the Yamada function built in Step 4. By the Ito formula, 

p,{sr) = P,(r) + 2 f p’,{sr)y^^d)l + f pUsy 2l'Lr.e ds + 2 f p'^iSl’^Sl’^ds. 

Jo Jo ^ Jo 

Proceeding as in Step 4, we find that 

( 33 ) pyin <^ + 2 £ py:nJydf3t + ^t + 


<r + 


C 


-t + 2jy'^iS:’nyy^dPs. 


l0g(l/£) JO 

Taking expectations, we deduce that E[(/Je(S'/’®)] < r + Ct/ \og(l/e), whence E[S'/’‘^] < r + £ + 
Ct/log(l/£). Coming back to (1551) and using the Doob inequality and that 0 <(/?/< 1, we 
conclude that 


E 


sup (/Je(S'[’^) 

[O.T] 


< r + 


CT 


l0g(l/£) 


+ C(rT + eT + 


T2 


log(l/£)) 


1/2 


< 


C(l + r)(r 


1 


log(l/£)) 


1/2 


because r G (0,1]. Then (I52ll follows from the fact that x < e + Pe{x). 
This allows us to conclude when y > Stt. 


□ 
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Proof of Theorem^^ when x > Stt. The existence of a (weak) solution to (1501) follows from Lemma 
m and the solution built there is frozen when it reaches 0. The pathwise uniqueness of such a 
frozen solution follows from the Lipschitz continuity of coefficients a, 6 on \ {0} and can easily be 
verified using the stopping times re = inf{t >0 : |^t| < 1/^} and that r = inf{t >0 : \Zt\ = 0} = 
sup^>]^ Ti (because 1 1 —>■ Zt is a.s. continuous on [0,oo)). Using Lemma [T51 we easily conclude that 
{Zf)t>o goes in law to this {Zt)t>o- Since and since the map z ha 

is continuous, we conclude that (Df)t>o goes in law, as e —>• 0, to (|Zt|“^/^Ztl{Zt/o})t>o- D 


To conclude the proof when y S (0, Stt), the only issue is to check the the uniqueness in law of 
the solution. We define h 2 Tr{ 0 ) = 9 — 2 ti[ 9 /G [0,27r). 

Lemma 19. Consider 0 < sq < to, a continuous function r : [0,oo) ha satisfying that 
= ^io = 0 rt > 0 for all t € (so,to) and J {rs)~^ds = oo for all t G (so,to)- There is 
a law r(so, ^0; (^s)sG[so.to]) C'(('S 0 ; ^o); [0, 27r)) (with the torus topology on [0, 27r)^ such that for 
any filtration (TLt)t>o in which we have a 1-dimensional Brownian motion (7t)t>o and a 

{'Ht)t>o-adapted process iTt)t(z(so,to) with Tt = /i27i-(Tlx + fjrs)~^'^^djs) for all sq < u < t < to, 
iTt)te{so,to) is independent ofUso T{sQ,tQ, {rs)s^iso,to])-distributed. 

Proof of LemmaHM Existence. Let uq G (soj^o) be chosen arbitrarily. We consider a Brownian 
motion independent of a random variable 0 , uniformly distributed on [ 0 , 27 r). We put Tt = 

h2Tv{Q + J^g(,rs)~^/‘^d'ys) for all t G (sofio) (with when t < uq). 

Then {Tt)t^(^so,to) is clearly continuous for the torus topology and it holds that Tt = ft. 27 r( 7 « + 
/^(rs)“^/^d7s) for all sq < m i < ^o- Furthermore, for each fixed t G (so,to), by independence 
between 0 and 7, the conditional law of Tt knowing ( 7 s)s>o is the uniform distribution on [ 0 , 27 r), 
which implies that Tt is independent of (7s)s>o- Finally, we have to verify that setting Ht = 
a{{Ts), (7s)se[o.t])! (7s)s>o is a (■Ht)t>o-Brownian motion. Let thus t G (soj^o) be fixed. We have 
to verify that (7s 7 i)sAt is independent of (T^, 7 s)se(so,i]* Since Tg — ^27r {Tt-Jg{ru) for 

all s G (soi^]: it holds that cr((T);, 7 s)se(so.i]) = a’fiTt, ( 7 s)sg(so.d) conclusion easily follows 

from the independence between Tt and (7s)se[so.to]■ 

Uniqueness. We thus consider a filtration {'Ht)t>o in which we have a Brownian motion (7t)t>o 
and an adapted process {Tt)t£(^so,to) satisfying T* = /i27r(T„ + /^(Ts)“fo^d7s) for all sq < m < t < to- 
We will show that for any fixed uq G (sq, to) I Tuq is uniformly distributed on [0, 27 r) and independent 
ofHso Vcr((7t)t>o). Since iTt)te[so,to] is cr(riio, (7* - 7sJtg(s(,_tj,))-measurable and since (7t)t>o is a 
(’Ht)t>o-Brownian motion, we conclude that iTt)t^(^so,to) is independent of TLso- Furthermore, the 
process (Tt)t^(^so,to) clearly has the same law as the one built above. 

For 0 < e < 7 ? < uo - So, we have T^^ = h27v{Ts„+g + X,*“(!7(Ts)“fo^d7s + Cg+r/(^s)~^'^^djs). By 
assumption, the vector (rs)~^''^d'ys, f^°^^(rs)~^''‘^d'ys) has independent components and is 

independent of Bsq V a(Tsa-i-e). Setting = fg°(!(g(rs)~^ds, we thus have, for any : R ha [ 0 , 00) 
continuous and 27 r-periodic, 


( 34 ) 


Ht{Tuo) I TLso V crLTso+e, (js - 

= [ ‘f(Tso+e+ [ {rs)~^^^d'ys 

«y M J Sq -j-T^ 


7soA'o)-s>soA'o)] 

\ e~^ 

+ as rix] — r^^dx —>■ ( 27 r)“ 



ip{x)dx 
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a.s. as £ —>• 0. This last convergence follows from the facts that lim£_>o = oo and that, setting 
if{x) := ip{x) - (27r)-i Jg"" (p{y)dy and ^{x) := /“ (p{y)dy, for all 6 G [0, 27r), 



^ x^[0,2'k) 


We used an integration by parts, that <i>(0) = $(271) = 0 and that |$(j/)| < 2?! sup 3 ,g[Q 2 ,r) 
for all 2 / G [0, 27r). 

We deduce from (El that Tuq is uniformly distributed on [0, 27r) and is independent of T-Lsg V 
cr(( 7 s — 'ysa+ri)s>sa+r])- Since r] > 0 can be chosen arbitrarily small, we conclude that is 
independent of TLso V cr(( 7 s - 'yso)s>so) = '^sa V cr((7s)s>o) as desired. □ 

Lemma 20. Assume that x G (0, Stt). There is uniqueness in law for (1301) among solutions sueh 
that a.s., l[Zt=o}dt = 0 . 

Proof. As in the proof of Theorem [T3 when x ^ Stt, (1301) admits a pathwise unique solution until 
it reaches 0. All the difficulty is thus to prove the uniqueness in law of the solution started at 
0 . We thus consider, if it exists, a continuous solution {Zt)t>o to (l30l) with Zq = 0 , adapted 
to some filtration {J-t)t>o in which (Wt)t>o is a 2 -dimensional Brownian motion, and such that 
l.^Zt=o}dt vanishes a.s. 

Step 1. We define Rt = |Ztp/^/4 and fit = l[Zs^o}\Zs\~^Zg ■ dWs, which is clearly a 1- 

dimensional (J^t)t>o-Brownian motion. Here we prove that 

(35) Ri = 2 f yJWsdPs + (2 - x/ (47r))t. 

Jo 

Starting from (j30p (with Zq = 0) and using the ltd formula, we easily find that 

\Ztf = 12 [ \Z,\^^^dfs + (72 - 3x/^) [ |Z«|4/3ds. 

Jo Jo 

For > 0, using again Ito’s formula, we find that 

{\Zt\^ + 72)1/3 =^1/3 + 4 r |Z,|3/3(|Z«|2 + y)-y^dfs + (24 - x/tt) T |Z,|4/3(|Zt|2 + 72)-2/3ds 
Jo Jo 

-16 [ \Zs\^°/^{\Zs\^ + v)-^/^ds. 

Jo 

Since /o ^{Zs=o}ds = 0 a.s. by assumption, the Lebesgue theorem ensures us that the sum of the 
two last terms in the right-hand side converges a.s. to (8 — x/^)/ as 7 —>• 0. The ltd isometry 
ensures that the second term in the right-hand side converges in to 4 f* \Zsf-/'^dfs. All in all, 
we find that \Ztfl‘^ = 4 \Zsf-/^dl3s -I- (8 — x/^)/- Dividing by 4 completes the proof of (1551) . 
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Step 2. We consider, for each 77 > 0, a nondecreasing C^-function : [0, 00 ) i-A [0, 00 ) such that 
■0(11) = 0 for all u G [0,77/2] and = 1 for all m > 77 . Observing that 

where '£',,( 2 :) = is of class on R^, we easily obtain, starting from (1511)) and 

applying the Ito formula. 


(36) 7/7^(i?7) 


A 

l^tl 


^P^{Rs){2Zi 


2Z., 


Z^.dW, 


|Z«|V3 |Z,|5/3 


*)+i 


\z. 


ii;'JRs)dRs + 2MRs)Rsds) 


where, for z G with respective coordinates zi and Z 2 , z-^ denotes the element of R^ with 
respective coordinates —Z 2 and zi. 

Let 7 t = /g l[z,^o}\Zs\~^Z^ ■ dWs- Since (/ 3 , 7 )t = fg lfz,^o}lZsl~^Zs ■ Z^ds = 0, the pro¬ 
cess ( 7 t)t>o is a 1 -dimensional (JT)t>o-Brownian independent of {Pt)t>o and thus also of (i?s)s>o 
(because (i?s)s>o is cr(i?o, (/3s)s>o)-nieasurable by pathwise uniqueness for the SDE it solves). 

For any 0 < m < t, on the event {inf[„ (] Rg > 0}, choosing 77 G (0, inf[^i_t] Rg) in the difference 
between (|36l) and the same equation with t replaced by u, we obtain 


(37) 


Z-t Zu f* Z^ d'js 2Zg \ Zu f* / Z^ d'jg Zg ds \ 


Step 3. For s > 0 such that i?s > 0 we define Tg G [0, 27r) through the equality \Zg\~^Zg = 

For s > 0 with Rg = 0, we simply put Tg = 0. We used the natural identification between R^ and 
C : for 0 G R, we denote by e*® (resp. ie*®) the 2-dimensional vector with coordinates coa0 and 
sin6) (resp. — sin0 and cos 9). We claim that for all 0 < it < t, on the event {inf[„ t] Rg > 0}, it 
holds that Tt = /72^(T„ -|- RJ^^'^d'ys). 

To check this claim, on the event {inf[„_(] i?s > 0}, we introduce =Tu + Rs for all 

V G [m, i\. Since (7«)«>o is independent of the event {inf[„ Rg > 0}, we can apply the Ito formula: 


for all V G [u,t], 



JTs 


ds \ 

^7 


Recalling (1571) and using a uniqueness argument, we deduce that on the event {inf[„ Rg > 0}, 
\Zy\~^Zy = whence Ty = h 2 TriTv) for all v G [u,t]. 

Step 4- Here we check that a.s., Rj^ds = 00 for alH > 0 such that Rt = 0 and all h > 0. 
This follows from the fact that for all T > 0, lim„ NO suptgjo j.] [u(lVlog(l/u))] 77? t+u 'n/ Rt I — 

•\/2 a.s., see Khoshnevisan [55J ( 2 . 1 a) p 1299] and recall that {Rs)s>o is a squared (2 — x/( 47 r))- 
dimensional Bessel process starting from 0 by Step 1, with 2 — x/( 47 r) > 0. 

Step 5. Here we verify that conditionally on (i?s)s>o, for any (T((i?s)s>o)“nieasurable finite 
family 0 < si < ti < S 2 < ^2 < ■ ■ • < Sn < such that for all k = 1,... ,77,, Rg,, = Rt^ = 0 
and i?s > 0 on {sk,tk), the variables {{Ts)se(sk,tk)^^ = l,...,n} are independent and for each 
/C = 1, ... , 77 , (Tg)gg (sfc.tfc) is r(sfe, tfc, (i?s)sG(sfc.tfc))“distributed. The function T was introduced in 
Lemma [121 

Let (Zt,g()t>o denote the canonical process on C'([0,oo),R^ x R) endowed with the conditional 
law of {Zt,jt)t>o knowing (i?t)t>o- We define T* G [0, 27r) by |Zt|“^Zt = e*'^* if Z* 7 -^ 0 and 
Tt = 0 else. We introduce the filtration Rt = tT((Tg, gg)sg[o_t])- We claim that a.s., (gt)t>o is a 
(’Ht)t>o-Brownian motion, because ( 7 t)t>o is independent of a{{Rg)g>o) and is a Brownian motion 
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in the filtration {J-t)t>o to which {Tt)t>o is adapted: for all t > 0, all bounded measurable <f>, ih, 
E ^'((7t+s-7t)s>o)^'((7s,T’s)se[o,t]) {Rs)s>o 
=E ^'(( 7 s,rs)sg[o,t])E ^>((74+8 - 7t)s>o) V cr((i?s)s>o) {Rs)s>0 
=E $((7t+s - 70s>o) E ^'((7s,Ts)sg[o,t]) {Rs)s>o ■ 


Fix now k € {1,... ,n}. It a.s. holds that Tt = h 2 ^(Tu + f* Rs ^^^dgs) for al\ Sk < u < t < tk 
by Step 3 and that Rj^ds = oo for all t £ {sk,tk) by Step 4. Applying Lemma fT^ we find 
that a.s., (Ts)sg(sfc,tfc) is independent of Rsk and is r(sfc, (i?s)j;g(j,^_(j,))-distributed. Using that 
is -measurable for each k = 1,... ,n, the independence easily follows. 

Step 6. By Step 1, {Rt)t>o is a (2 — x/(47r))-dimensional Bessel process starting from 0. By 
Step 5, the conditional law of knowing {Rt)t>o is also determined: conditionally 

on (i?s)s> 0 ; for any (T((i?s)s>o)-measurable finite family {{sk,tk),k = l,...,n} of excursions of 
{Rs)s>o, we know the law of iTs)seu^^^isk,tk)- Since by construction Zt = (4i?t)3/2e*^‘ 
the law of {Zt)t>o is thus entirely characterized. □ 


Finally, we can give the 

Proof of Theorem^T^ when x £ (0, Stt). First, the existence of a solution {Zt)t>o to (l30l) such that 
a.s. l[Zt=o}dt = 0 follows from Lemma [TSl the solution {Zt)t>o built there satisfies that 

|^^|2/3/4 is a (2 — x/(47r))-dimensional Bessel process, whence l[Zt=o}dt = 0 a.s. by [35J p 
442]. The uniqueness in law of this solution has been checked in Lemma [201 The convergence of 
{Zf)t>o to {Zt)t>o clearly follows from Lemma |18l and from this uniqueness in law. This implies 
as in the case x > Stt that {D^)t>o goes in law to {\Zt\~^^^Ztl^Zt=/=o})t>o- 

It remains to verify that when x £ (0,47r), Dt = solves (^5)) and that unique¬ 

ness in law holds true for (l28l) . 

For {Dt)t>o a solution to (l28l) . one easily checks by Ito’s formula that Zt = \Dt\‘^Dt solves (l30l) 
and that \Dt\^ is a (2 — x/(47r))-dimensional Bessel process, whence l^Zt=o}dt = 0 a.s. by [321 
p 442]. The uniqueness in law for (1331) then follows from Lemma 1301 

For {Zt)t>o built above, by Ito’s formula, for rj > 0, 

(|Z*|2 + r,)-^/^Zt = i\Zo\^ + + 2 [ |Z,|2/3(|Z,|2 + rjY^/YWs 

Jo 

+ 4^ (\ZsY^^H\Zs\^ +vY^^'' -\Zs\^^Y\Zs\^ + r]Y^^^)Zsdl3s 
+ ^ (^(16 - 3x/(27r))|Z,|-2/3(|Zt|2 + r,Y^/^ + (xA - 48)\Z,YY\ZY + vY^^" 

+ 32\ZsY/Y\zY + pY^^^^Zsds. 

By the Ito isometry and the Lebesgue theorem and since a.s. = 0i the second term 

on the RHS tends to 2IUi in and the third term on the RHS tends to 0 in L^. Since |Ztp/^/4 
is a (2 — x/(47r))-dimensional squared Bessel process and 2 — x/(47r) > 1, [321 Exercise 1.26 p 451] 
ensures that a.s. < oo. Hence the Lebesgue theorem ensures us that the last term on 
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the RHS converges a.s. to —(x/(27r)) \Zs\ We conclude that Dt = \Zt\ Ztl^Zt=^o} 

solves Dt = Do + 2Wt — {x/{2n)) \Ds\~'^Dsds, which completes the proof. □ 

7. On the system with N >3 particles 

7.1. Classification of refiecting and sticky collisions. We have seen in the proof of Lemma 
Ilhl Step 2 that very roughly, the empirical variance of the positions of k particles in the system 
with N particles resembles a squared Bessel process of dimension Sn,x{^) = 1)(2 —X^/(47rA^)). 

Fix X > 0 Slid N > 3 and consider the regularized particle system which is always well-posed. 
We now describe formally the expected behavior of its limit as e —>• 0. According to [3U Page 442] 
and the comparison theorem [32l Theorem 3.7 p 394], the following events should occur: 

• if > 2, no collisions of subsystems of k particles, 

• if 6N,x{k) S (0,2), (instantaneously) reflecting collisions of subsystems of k particles, 

• if 6N,x{k) < 0, sticky collisions of subsystems of k particles. 

Let us now study the inequality ^jv,x(fc) ^ 2. We have already seen in the proof of Lemma 
[Ej-Step 2 that when x G (0, 87r(A^ — 2)/{N — 1)], (5jv.x(^) > 2 for all k G {3,...,A^}. When 
X G — 2)/{N — l),47rAf/3], (57v,x(3) > 2 whereas (5jv.x(2) < 2 and 5n,x{N) < 2, hence 

the two roots ^ = [1 -I- (87rA^)/x ± y^(l -I- 8ttN/xY — 647rA^/x]/2 of the second order equation 
5n,x{x) = 2 are such that xjj ^ G (2,3] and x'lf ^ G [3,Af), so that 5n,x{‘^) < 2, 5n,x{^) — 2 for 
fc G {3,..., [x~l^ } and Sjx^xi^) < 2 for fc G {+ 1,..., A^}. Finally, one easily checks that 

xjj — d and iT^N/s ~ d- By strict monotonicity of the map x ^N,x{k)j Y^e conclude that 
if X > 47rA^/3, then SN,x{k) < 2 for all fc G {2,..., N}. 

Let us next study the inequality SN,x{k) < 0, which, for fc G {2,... ,N} is equivalent to fc > 
SttN/x- Hence for x G (0, 87r), SN^^{k) > 0 for all k G {2,...,N} whereas for x G [87r,47rA^), 
dN,xik) > 0 for all /c G {2,..., \8ttN/x] — 1} and dN,xik) < 0 for all k G {[STrAf/xl,.. ■, N} with 
the two sets non empty. When x > 4:TtN, SN,x{k) < 0 for all fc G {2,..., N}. 

When N > 6, we end up with the following picture. 

(a) If X £ (0, 87r(A^ — 2)/{N — 1)], the regularized particle system should tend to the particle 
system O and the latter should have a unique (in law) solution. Indeed, it holds that 6N-,x{k) > 2 
for all fc > 3 and that 5n,x{‘^) G (0,2), so that only binary reflecting collisions occur. We have 
already checked a tightness/consistency result in this spirit in Theorem [7| Only the uniqueness in 
law remains open. 

(b) If X £ {8tt{N — 2)/{N — I),87r), the regularized particle system should tend to the particle 

system O and the latter should also have a unique (in law) solution. One may check that fco •= 
L^ivxJ + 1 S ~ (il suffices to verify that 5n,x{^ — 2) > 5n,s,-k{N — 2) > 2). In this 

situation, there should be binary reflecting collisions and also reflecting collisions of subsystems of 
particles with cardinality in {fco, A^}- To check the existence (and a fortiori uniqueness) of such 
a process, one has to control the drift term during the collisions with reflection. In the present 
paper, we are more or less able to contol the drift during a (reflecting) binary collision, but we 
have not the least idea of what to do during a fc-ary refiecting collision with k > 3. 

(c) If X G [Stt, 47rA^/3], the regularized particle system should tend to a particle system with 

sticky collisions that we will describe more precisely in the next subsection. One can check that, 
for ko := -|- 1 > 4 and fci := [STrAf/xl < N, we have fco S {fci — 2, fci — 1} (just verify that 

SN,x{ki — 3) > Sn^^{8ttN/x — 2) > 2). Thus, binary reflecting collisions, as well as fc-ary reflecting 
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collisions, for k G {kg, ki — 1}, should occur, as well as sticky collisions of subsytem of A:-particles, 
for k G {ki,..., N}. Assume e.g. that ko = ki — 1. What might happen is that, at some time, kg 
particles become close to each other, they may collide (with reflection) a few times, then another 
particle is attracted in the zone, the ko + 1 = ki particles meet and then remain stuck forever. 
Such a cluster will move with a very small diffusion coefficient and should collide later with other 
particles (or clusters) in a sticky way. Of course, such a result would be very interesting but it 
seems very difficult to prove, because to check the existence of such a process, one would have to 
control the drift term during the collisions with reflection, as mentioned previously. The sticky 
collisions should be easier to describe. 

(d) If X € (47rA^/3,47rA^), the same situation as previously should arise, except that there 
should be fc-ary reflecting collisions for all /c € {2,..., [STriV/x] — 1} and sticky fc-ary collisions for 

allfcG{r8^iV/xl,...,iV}. 

(e) If finally x > 47rA^, then there should be sticky k-ary collisions for all k G {2,..., iV}. 

When N = 5, we find the following dichotomy. If x S (0, 67r], only binary reflecting collisions. If 
X G (Gtt, 207r/3], only reflecting collisions of subsystems oi k G {2,5} particles. If x £ (207r/3, Stt), 
only reflecting collisions of subsystems oi k G (2,3,4, 5} particles. If x £ [Stt, 207r), reflecting 
collisions of subsystems of /c G (2,..., |’407r/xl — 1} particles and sticky collisions of subsystems of 
k G { [dOTr/xl,..., 5} particles. If x ^ ^Ott, fc-ary sticky collisions for all G (2,..., 5}. 

When finally N G {3,4}, 8'!t{N — 2)/{N — 1) = AttN/H and the situation is as follows. If 
X G (0, 8Tr{N — 2)/{N — I)], only binary reflecting collisions. If x S (A-kN/Z, Stt), reflecting collisions 
of subsystems of fc G {2,...,A^} particles. If x t: [87r,47rA^), fe-ary reflecting collisions for k G 
{2,..., |'87rA^/xl — 1} and fc-ary sticky collisions of subsystems of fc G {{87riV/xl, ■. ■, JV} particles. 
If X > 471A^, fc-ary sticky collisions for all fc G {2,..., iV}. 


7.2. A particle system in the supercritical case. When x > Stt, the following dynamics 
should describe the limit of the regularized particle system as e —>■ 0. Particles are characterized 
by their masses and their positions. Initially, we start with A particles with masses Uq, ..., all 
equal to 1/A and with some given positions Aq’ ,..., Ag ’ . If now at some time t > 0, we have 
At particles (At will be a.s. nonincreasing) with masses (such that = l)i '^6 

make the positions evolve according to 


(38) 


dXl’^ = 


Nvl 


dBl 


Nt 


?A(A, 


i,N 


-Xi 


,N 


)dt, 


= 1 , 


,At 


until the next collision between at least two of these At particles. If the sum S of the masses 
of the particles involved in the collision is smaller than 87r/x, they should automatically separate 
instantaneously and we carry on making evolve the system according to (I38|) (with the same values 
for the masses and for At) until the next collision. If now S exceeds 87r/x, the particles involved 
in the collision are replaced by a single particle with mass S, the number of particles is decreased 
accordingly, the particles are relabeled, and we make evolve the system according to (l38l) with 
these new values for At and for the masses until the next collision. 


By construction, the masses take values in {I/A, 2/A,..., A/A} and actually in {fc/A : fc = 1 
or SttN/x < k < N}. A particle of mass fc/A with fc > 2 has to be seen as a duster oi k elementary 
particles. The drift term is thus easily understood: a single elementary particle interacts with the 
other ones proportionaly to 1/A, so that a cluster consisting of fc elementary particles interacts 
with the other ones proportionally to its mass fc/A. The diffusion coefficients are also quite natural: 
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a single particle being subjected to a Brownian excitation with coefficient a cluster with mass 
k/N is excited by the mean of k Brownian motions with coefficient y/2, that is, by a Brownian 
motion with coefficient ^/2jk. 

If Nt > 2, setting for I C iVt} with cardinality |/| > 2, X/ = 

and Rf = iN/2)J2i^i ~ a simple computation shows that, when neglecting the 

interaction with particles with label outside I, behaves like a squared Bessel process of dimension 
2(|/| — 1) — jei i^j ^ (l-^l ~ 1)[2 — ‘5'x/(47r)], which is nonpositive as soon as 

S = ^ Stt/x- 

Let us mention that once a cluster is formed, its mass necessarily exceeds Stt/x, so that any 
collision involving a cluster will be sticky. 

The existence of such a process is not clear. Sticky collisions should not be very hard to treat. 
The main difficulty is to control reflecting collisions. As explained just above, reflecting collisions 
only concern particles with masses 1/N, so that the classification given in Subsection 17.11 should 
still be relevant. Thus we believe that the main difficulty is to build a (necessarily nontrivial) local 
(in time) solution to (|T|) when x ^ Stt and starting from an initial condition where k particles have 
the same initial positions, for some k G {2,..., \8^^N/x^ — !}• 

7.3. Comments. Observe that this process is different of the one introduced by Haskovec and 
Schmeiser in [10] where they consider a system of particles with different masses to approximate 
the singular solution to the Keller-Segel equation. In fact, rather than considering like us the 
limit e —>■ 0 of the regularized particle system ([9]), they first prove propagation of chaos as X —>■ 
oo for a fixed £ > 0 in m- More precisely, they check that for fixed A: > 1, the density of 
,..., solving ([9]) (with another regularized kernel K^) converges as X —>• oo to 

rii^i /f where (/f)t>o solves the regularized Keller-Segel partial differential equation 

dtftix) + x<i^^xiiKe * ft){x)f^ (x)) = 

The limiting behaviour of (/f )t>o as e —>■ 0 was studied in [4] and involves a defect measure. Then 
Haskovec and Schmeiser introduce in |10j a particle system associated with this limit, in which 
there are heavy particles that occupy a positive proportion of the mass, interact with the other 
particles, but do not undergo any Brownian excitation. 
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